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1. Introduction 

Recent work on unstable DO-branes of two-dimensional bosonic string theory [jl|-|^ has 
led to reinterpretation of the well-known large- matrix quantum mechanics formulation 
of this theory (for reviews, see ||^-|T^) as exact open/closed string duality. The open 
strings live on unstable DO-branes; the boundary state of such a DO-brane is a product 



of the ZZ boundary state [|T3[ for the Liouville field, localized at large (j), and of the 
Neumann boundary state for the time coordinate [^,^ . The dynamics of these open strings 
is governed by a gauged quantum mechanics of an x Hermitian matrix M with an 
asymmetric (e.g. cubic) potential. This model is exactly solvable since the eigenvalues act 
as free fermions. In fact, these eigenvalues are the DO-branes. In the double scaling limit 



1^ , the ground state of the 2-d bosonic string theory is constructed by filling one side 
of the inverted harmonic oscillator potential, — A^/(2q;'), with free fermions up to Fermi 
level —11 as measured from the top of the potential. Since Qs ~ this state has a 

non-pertubative tunnelling instability. 

While this matrix model formulation of the 2-d closed bosonic string has been known 



for quite some time fT^, similar formulations of NSR strings have been a long-standing 
problem. In recent work [lT6| , p!7[] a solution of this problem was found for two-dimensional 
type strings. Let us briefly summarize the logic that led to this solution. Consideration 
of unstable DO-branes of the type OB theory indicates that the dynamics of open strings 
living on them is again governed by a gauged Hermitian matrix model, but now with a 
symmetric double- well potential. This led the authors of [IB, 171 to conjecture that the 
ground state of the closed OB string theory corresponds to filling the potential — A^/(4a') 
symmetrically up to Fermi level — /U. In the continuum formulation the parameter ^ enters 
the superpotential of the super-Liouville theory as ^ie'^ . This explains why in the OB theory, 
unlike in the bosonic string, can have either sign. In fact, this theory has a symmetry 
under n —fx which was called S-duality in For either sign of ^ the fermions 

are divided symmetrically into two branches in phase space; hence, we may loosely call 
this a two-cut eigenvalue distribution. 

A matrix model formulation of the 2-d type OA closed string may be derived in a 
similar fashion. The OA theory has charged DO-branes and anti DO-branes. The dynamics 
of open strings on N + q DO-branes and A^ anti DO-branes is described by (A^ -t- q') x A^ 
complex matrix quantum mechanics with U{N + q) x U{N) gauge symmetry. Just as 
the OB model, this model is exactly solvable in terms of free fermions and is stable non- 
perturbatively [|17 . 
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In [16-23] these models were studied further. In particular, the matrix model formu- 
lations of the 2-d type strings were subjected to a number of stringent checks vs. the 
continuum worldsheet formulation in terms of c = 1 super-conformal field theory coupled 
to super-Liouville theory. In this paper we consider further extensions of these dualities to 
c < 1 theories coupled to super-Liouville. If we turn on relevant operators in the c = 1 the- 
ory, approporiately dressed by the Liouville field, then the theory undergoes gravitational 
RG fiow to c < 1 models coupled to the super-Liouville theory. Therefore, we expect such 
string theories to have matrix model duals closely related to the ones found for c = 1. In 
this paper we indeed argue that type OB theories are again dual to double-scaling limits 
of hermitian matrix models with two-cut eigenvalue distributions (or, equivalently, of the 
unitary matrix models), while type OA theories are dual to complex matrix models. 

For simplicity, we will restrict our discussion to the one-matrix models. Large 
unitary matrix models of this type were solved in ||2^ - |30|] while the complex matrix models 
31-35|. The generic critical behavior of such models is that of the c = theory (pure 



m 



supergravity) , i.e. of type strings in one dimension. The unitary matrix model below 
the Gross- Witten phase transition has a two-cut eigenvalue distribution and above the 
transition it has one cutJl In the double scaling limit we conjecture that this model is dual 
to OB string theory with n > and ^ < respectively. The complex matrix model has a 
phase transition where the eigenvalues reach the origin. We conjecture its double scaling 
limit around this phase transition to be the dual of OA string theory. It was observed long 
ago that the double scaling limit of the generic complex matrix model is equivalent to the 
one of the generic unitary matrix model Indeed, for c = we find that a change in 
the sign of the left-moving fermion on the wordsheet sends the OA theory at some value of 
H to the OB theory at — /i. 

By fine-tuning the potential in the matrix integrals, we are able to also describe a 
certain class of non-unitary c < SCFT's coupled to super-Liouville theory. A worldsheet 
interpretation of these matrix models has been a longstanding puzzle. In fact, the idea 
that the two-cut hermitian matrix models are dual to SCFT's coupled to supergravity was 
advanced in the early 90's |l36|j27| , p8| but, as far as we know, was not tested thoroughly. 
We will present a number of consistency checks of this duality conjecture which rely on 
the interplay of the OB and OA models. 



^ Later in the paper we will sometimes use the phrase "two-cut hermitian matrix model" in 
referring to the model around this transition, either above or below the transition where it has 
one or two cuts. 
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In section 2 we discuss the unitary matrix models, reviewing and extending the existing 
hterature on the subject. Section 3 is devoted to various aspects of the matrix model 
resolvent. It satisfies a loop equation which describes a Riemann surface. We demonstrate 
the discussion with the explicit solution of the simplest model at tree level. The Riemann 
surface also leads to a new insight into the FZZT []37|j3^] and the ZZ [|T3| branes in the 
theory. In sections 4 and 5 we present a detailed analysis of the simplest nontrivial theory 
- c = theory - pure supergravity, and the first multicritical point. In section 6 we turn to 
discussion of the complex matrix model and its solutions. Section 7 presents a "spacetime" 
picture of these models. In section 8-10 we discuss the worldsheet interpretation of these 
theories. We analyze the R-R vertex operators and uncover interesting dependence on 
the sign of the cosmological constant, we discuss the properties of superminimal models 
coupled to supergravity and we explore the torus amplitude in these theories. In section 
11 we present our conclusions and open questions for future research. 

Several appendices provide more details for the interested reader. Some of these 
details are reviews of known results. In appendix A we discuss the comparison between 
the worldsheet and the matrix model results for the first multicritical point. In appendix 
B we mention some properties of superconformal minimal models and a restriction due to 
modular invariance on such models. Appendix C includes an assortment of results about 
the Zakharov-Shabat hierarchy of differential operators. In appendix D and appendix E 
we discuss various aspects of the complex matrix model. 



2. Unitary Matrix Models 

In this section we study unitary matrix models. The unitary one-matrix integrals have 
the form 

Z = j dUexp(^-^TtV{U + U^)^ , (2.1) 

where t/ is a unitary Nx N matrix. The simplest such model, with potential ~ Tr{U+U^), 
is obtained if one considers Wilson's lattice action for a single plaquette. This one-plaquette 
model was originally solved in the large- limit by Gross and Witten . For 7 = oo the 
eigenvalues e*^ are uniformly distributed on the circle parameterized by 9. As 7 decreases, 
the eigenvalue distribution gets distorted: it starts decreasing in the region where the 
potential has a maximum. As 7 is decreased below a critical value, 7c, a gap appears in 
the eigenvalue distribution: this is the third-order large phase transition discovered in 
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2|]. The generic (and simplest) case is when the potential has a quadratic maximum. 



This is labelled as /c = 1 in the classification of critical points. By further fine-tuning 



the potential, Periwal and Shevitz ||2^ found an infinite sequence of multi-critical points 
labelled by a positive integer k, and found their descriptions in the double-scaling limit in 
terms of the mKdV hierarchy of differential equations. 

The double-scaling limit of the simplest unitary matrix model, k = 1, with potential 
Tr{U + U^) is described by the Painleve II equation ||2^ 



2r-f + xf = 0, (2.2) 

where x ~ (7c - 7)A^^/^. The free energy F{x) is determined by F" = /(x)^/4. We will 
identify this model with pure 2-d supergravity where x is proportional to the parameter /x 
in the super-Liouville interaction, and —F{x) is the sum over surfaces. 

More generally, the k-th critical point is described by a non-linear differential equation 
for f{x) of order 2k. The solution / has the following large x expansion: 

fix) = xi/^^'^) (^1 - ^^^.-C^'^+D/'^ + , (2.3) 
from which it follows that 

^{2k+i}/k 2k + 1 

= - 4(2 + lA)(l + lA) - ^ '"^ + • • ■ ■ '2-^' 

In the double scaling limit the physics comes from the coalescence of two cuts. So, 
a double scaling limit of a hermitian matrix model around the point when two cuts meet 
will lead to the same free energy as the unitary matrix model. Important further steps 
in studying these models were made in p^ , p8| , [29[ 1 where both odd and even perturbations 
to the potential were considered. The complete model is described by two functions, f{x) 
and g{x), which in general satisfy coupled equations. An integer m specifies the critical 



points studied in (for even m the relation to k of is m = 2k). In the m = 2 case 
the equations are 

2f" + f{g'-f)+xf = 0, 

(2.5) 

2g" + gig^-f) + xg = 0, 
while the free energy is determined by 

F" = \[f{x)'-g{x)'] . (2.6) 



More generally we will need two functions F^, Gi which are polynomials in f, g and 
their derivatives. These polynomials are related to the Zakharov-Shabat hierarchy, which 
generalizes the mKdV hierarchy. They are defined by expanding the matrix resolvent 



^ oo 

O = {x\J^R\x) = {x\———-—\x) = V (-JiF, - iJ2Gi + JsHi)C^'+^^ (2.7) 

U + - C,J3 

where Ji = ai/2 with ai are the standard Pauli matrices and 

Our problem is invariant under a "boost" symmetry of / and g. Therefore, it is 
natural to define 



/±^ = -"^ Q= 'l']=r[:-s . (2.9) 



f-g ; "^[e-^ 
Similarly, we define a new "boosted" operator 

oo 

d = e-^^^Oe^-^^ = J2 i-JiRi - iJ2^i + J3Hi)C-'^^+^^ (2.10) 
i=-i 



where Ri and are related to Fi and Gi by a boost. Using ( |2.1C1| ) and ( |2.7D we find 



Q= 7 If 7TT' ^^^^ '^ = ^'- ^2.11) 

D + rJi + [lv - QJs 

We see that undifferentiated (3 does not appear in O, only u and its derivatives appear. 
We can see from that a constant shift in w results in a shift in ^, which in turn 

produces a redefinition of the expansion where each Hi gets mixed with lower Hi terms. 
The 2x2 matrix O obeys the equation 

[O,D + Q-(:Js]=0 (2.12) 

which leads to recursion relations for F;, Gi and Hi in ( pj.T] ). So we can determine them all 
from the lowest ones G-i = F-i = 0, H-i = 1 [^. Equivalently, we can derive recursion 
relations for Ri, G^. It turns out that rQi is always a total derivative. In fact, rG; = —H^. 
The recursion relations then become 

Ri+i = uRi - (^] + rHi 



\ r J (2.13) 
= uHl - rR[ . 
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The first few terms are 



= 1 , = ; 

ifo = , Ro = r; 
Hi = , Ri = ujr ; 

H2 = -r^uj , R2 = + ruo^ + r" ; ^^'^^^ 

84 322,1/2 // 

H3 = —r r uj -\ — r — rr , 

3 

R3 — —-r^u + ru^ + Sr'uj' + "iujr" + ruj" . 
These equations may also be obtained by "boosting" the recursion relations for Fi, Gi in 



It is an interesting exercise to find the terms in Hi and Ri with no derivatives of r or 
cj-ll In terms of the variables 



p2=^2_^^2^ COS99=^, (2.15) 



we show in appendix C that 



2 



Hi = ~ p^+Mcosv?Pi(cosv?) - Pz+i(cosv?)] = -p^+^^-^-^Pl{cos^) , 



The string equations of the matrix model may be stated as 



+ ^)tiGi = = + (2.17) 



l>0 1>Q 



where to ~ We also assume that = for / > m. These equations follow from varying 
the actioni 

/m 
dxY,tiHi+i. (2.18) 
1=0 

The l-th term correspond to a perturbation in the potential of the schematic form V ~ 
— tiZ^Tr[M^+^]. As pointed out in ||2^, the terms with odd / (and real ti) correspond to 
imaginary terms in the potential. Since we are doing an integral, there is nothing wrong 



^ Since uj = /3' this includes terms with derivatives in terms of tiie original variables. 
^ This can be proven as in the KdV case. See |1C] for a nice discussion of the KdV case. 
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with having these imaginary terms. In fact, they make the matrix integral more convergent, 
otherwise an odd term in the potential would be unbounded below if it ever dominates.S 
In all these models the free energy obeys F" = —Hi/2 = The action ( p.l8|) is 

invariant under x-independent shifts of /3. Therefore, the equation of motion of /3 is a total 
derivative and can be integrated by adding an integration constant q. Alternatively, we 
can view u — (3' as the independent variable in ( p.l8|) , and add to it qu: 

/m 
dx{J2tiHi+i+quj). (2.19) 

1=0 

If we assume that this action can be used in a quantum theory such that is well 
defined, invariance under (3 ^ (3 + 2m leads to the conclusion that q must be quantized. 
In our case we are only solving the classical equations coming from this action, and it is 
not clear to us which problem the quantum theory is the answer for. Below we will discuss 
the physical interpretation of q. 

The action ( |2.19| ) leads to the equations of motion 



S 

du{x 



/m m 
dxC^tiHi+i + quj) = -^ti{l + l)Ri = , 

1=0 1=0 

/m m 
dxC^tiHi+i+quj) = ^ti{l + l)Hi + q = Q . 
1 r\ 1 n 



(2.20) 



1=0 1=0 



Suppose we consider the m-th model, which has fixed t^. At first it seems that there 
are m parameters that we can vary: ti with / = 0, . . . m — 1. However the ability to shift 
w by a constant shows that we can set tm-i = at the expense of an analytic change of 
variables for the rest of the ti\ i.e. the operator that couples to t^n — 1 is redundant (this 
is analogous to a similar operator in the bosonic string, which was discussed in [^). As 
a result, there are only m — 1 operators. If we assign dimension minus one to x ~ to? 
then the operators have dimensions 1 — l/m for / = 0, ■ ■ ■ , m — 2. We will later match 
these operator dimensions with gravitational dimensions in (2, 2m) super-minimal models 
coupled to super-Liouville theory. 

If we set ti = for all odd / then there exists a solution with g = 0. We find that 
F2k{f,g = 0) is the k-th member of the mKdV hierarchy derived for the unitary matrix 



If we made tm imaginary for odd m we would have a real potential. Then the string equation 
becomes real if we define g ^ ig. But in this case there is no real solution for /, g 
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models in [25]. This corresponds to having an even potential and considering only even 
perturbations. For these even models m = 2k. We present examples of models with even 
m in sections 4 and 5. 

The simplest example of an odd model is the m — 1 theory. In this case ( p.l9|) becomes 

/I I 
dx{--r^uj - -xr^ + quj). (2.21) 
z z 



The equations of motion of this theory are 

{u + x)r = 
= 2q 



(2.22) 



and therefore u = = q/2. This "topological point" has been discussed in and 



we now interpret it as associated with nonzero q. We will say more about the physical 
interpretation of q below. Other examples of odd models are discussed in Appendix C. 



3. Comments on the Matrix Model Resolvent 

In this section we will analyze the matrix model in the standard large 't Hooft 
limit, i.e. in the planar approximation. Our main goals will be to identify the meaning of 
the parameter q introduced above, as well as understand the relation between the FZZT 
3^,^ and ZZ Jl^ branes in string theory. 



Consider the hermitian matrix model with potential V{M). We will be interested in 
two closely related operators. The macroscopic loop (the FZZT brane) 

W{z) = - — Trlog(M -z) = lim U j^Tre^i^'^^ + logej (3.1) 

(e is a UV cutoff) and the resolvent 

, dW(z) 1^1 , , 

^ W = ^ = TrTr— . 3.2 

Without causing confusion we will denote by W and R both the matrix model operators 
and their expectation values in the large theory. We will later think of R{z)dz = dW{z) 
as a one-form. It is clear from the expression for W{z) that it can have an additive 
ambiguity of with integer k. If this is the case, the one-form R{z)dz is not exact. 



Using the invariance of the matrix integral under SM = we derive the loop 
equation 

1 TV 1 U^/^miUo (3.3) 

^ M - z M - z I \ M -z I ^ ' 



The last term can be written as NV'i^z) (Ttj^^J plus a polynomial of degree n — 2 if 
the degree of V is n. In the large N limit the first term factorizes and we find the loop 
equation for the resolvent 

R{zf + r{z)R{z)-^f{z)=0, (3.4) 

where f{z) is a polynomial of degree n — 2. The solution of ( |3.4| ) is 

2R{z) = -V'{z) ± y/V'izy + f{z). (3.5) 

The cuts in this expression mean that it is a function on a Riemann surface which is a 
two-fold cover of the complex plane. For each value of the parameter z there are two points 
on the Riemann surface 

y' = V'{z)' + f{z). (3.6) 

We will denote them by P±{z). They differ in the sign of y, y{P^{z)) = —y{P-{z)) = 
^JV'{zY' + f{z)- The asymptotic behavior as z — >• oo should be i?(P+(z)) ^ -i) in the 
upper sheet. This determines 

2R{P±[z)) = ±^/V'{zr + f{z) - V'{z). (3.7) 

In the planar limit the eigenvalues form cuts. Their density p{\) = jj J2i <^(A — A^) is 
supported only on the cuts and it can be used to compute expectation values of operators, 
e.g. 

R(z) = — / Tr — - — ) = f dX-^^ for z not on the cuts (3.8) 

N \ M — z I J A — z 

from which we find 

Rix + ie) + Rix - ie) = 2V [ dX-^^ for x on the cuts, (3.9) 

J X- X 

where V J denotes the principal part. 

Using the expression for the integral in the planar limit we define the effective potential 
of a probe eigenvalue at z away from the cuts 

Veffiz) = V{z) -2 J dXp{X) log{z - A) = V{z) + 2W{z). (3.10) 



The force on a probe eigenvalue away from the cuts is 



Feffiz) = -V^ffiz) = -V'{z) - 2W'{z) = -V'{z) - 2R{z) 

(3.11) 

= -R{P+{z))+R{P.{z)) = -y{z). 

The force on an eigenvalue on the cut, — | {^effi^ + + ^effi^ ~ ^^)) vanishes, as can 
be verified using ( |3^ ) and 



V'{x) + 2V [ dXj^ = V'{x) + R{x + ie) + R[x - it) = 0. (3.12) 



We can also analytically continue the FZZT brane ( |3.1j ) to the second sheet and dis- 
tinguish between W{P^{z)) and W{P-{z)). Since the analytic continuation of an analytic 
function is unique, this generalizes the analytic continuation discussed in [^0]. 

This discussion makes it clear that the force ~y{z)dz or Rdz are one forms on the 
Riemann surface, while W is the corresponding potential. We can calculate the periods of 
the force around the a-cycles of the Riemann surface 

— I Rdz = -— (3.13) 

where Ci is a contour around the cut i and Ni is the number of eigenvalues in that cut. If 
there is only one cut, A'^i = A^. More generally, the residue of i? at z ^ oo is — 1. This 
discussion makes it clear that W{z) of ( |3.1|) indeed has additive ambiguities of 2'Kik/N 
with integer k. 

Similarly, we can study the periods of Rdz around the 6-cycles of the Riemann surface 
( p.6|) . Assume for simplicity that there is only one such cycle; i.e. there are only two cuts. 
It is a contour which connects the two cuts in the upper sheet and then connects them 
also in the lower sheet. The period is 

^ - ^ ^ydz = ^ (p Rdz. (3.14) 



47rz Jf, 2ni 

This means that W has an additive monodromy of 2niq/N. We will later identify it with 
the flux or the number of D-branes. This discussion suggest that q must be quantized. 
It does not prove it because this discussion applies only in the large limit with flnite 
q = q/N which is arbitrary. 
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3.1. Solution of the simplest model in the planar limit 

Now we will consider the simplest double- well model with V'{z) = ■^{z^ — z) which 
leads to the k = 1 critical behavior. The second order polynomial f{z) is determined by 
two parameters: 

f{z) = --z^ + hz + fo , (3.15) 
9 

where the first term is fixed by the condition that R{z) —>■ —1/z for large z. We impose 
that the sixth order polynomial y'^{z) has one double zero at —iuj {uj can in principle be 
any complex number; we will later take it to to be real) and four simple zeros. This gives 
one relation between the remaing two parameters and leads to the curve 

y'^ = l^(z + iujf iz^ - 2iuj z^ -z^(2 + 3lj^) +4iLj z(l + uj^) + A) , (3.16) 

with 

A = (1 - Ag + Gu'^ - 5uj'^) (3.17) 

Turning on the remaining parameter would spread the zero of y into a branch cut. 

The scaling limit is obtained when two of the simple zeros of y'^ go to infinity. Equiv- 
alently, we may scale z and u) to zero at fixed ratio, while neglecting higher powers of z 
and LO. It is clear from the curve ( |3.16| ) that A should scale like z'^. Therefore, the critical 
limit of g is g = 1/4. So, we parametrize g = | — x/2, with x scaling as u. In this limit 
the curve (|3.16|) becomes 



gly^ ^ -2{z + iuf {{z - iuf - 4u) (3.18) 

with gc = 1/4 and 

4u = x + 2uj^ (3.19) 

It is straightforward to compute the period of this y around the cut between the simple 
zeros z± = iiv ± I^Ju 

^^jydz = ^V2uu = q= ^. (3.20) 

It turns out that the second derivative of the free energy is d^F = 8u. Up to a rescaling 
of u, UJ and x, these are the same as the equations that follow from ( |2.2C1| ) for the simplest 
case, t2 = 1/3, to = 3:/2 and q nonzeroi. We will discuss these equations in more detail in 
section 4. 



^ The rescaling is Xthere = 2Xhere, Uthere = 2Uhere and Wthere = V^Where- 
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Now we solve the two equations (|3.19|) ( |3.20|) . One can identify u with the second 
derivative of the free energy. We wiU discuss first the case of q = and then q ^ 0. For 
g = 0, either u = or uj = 0. With a; = ( |3.16|) has a double zero at the origin and two 
simple zeros at z± = ±2y/u. It is natural to run the cuts in R from these two zeros to 
infinity. We take x = 4u to be positive and identify this with the positive x phase of the 
two cut model. The zero of y at the origin shows that an eigenvalue at that position feels 
no force. In string theory such an eigenvalue corresponds to an unstable ZZ brane. The 
potential has a maximum at this point. 

With u = (|3.16|) has double zeros at ±iuj. Here we take u to be real to describe 
the negative x = —luP' phase of the two/one cut model. The ends of the other cut of this 
model were scaled to infinity and are not visible in this limit. The effective potential does 
not have stationary points for real zq. Therefore it is not clear whether the system has 
unstable ZZ branes. However, one can place some eigenvalues at ±zc<; where the effective 
potential is stationary. We will see that this is equivalent to taking q ^ ^ and will be 
interpreted in the string theory as a charged ZZ brane. As we will soon see, branes at 
zw, along the positive imaginary axis, are charged D-branes and branes at — zw, along the 
negative imaginary axis, are charged anti-D-branes. 

Now let us turn to the case that q' 7^ 0. We take q — > 00 with finite g = so that its 
effects are visible in the planar limit. 

Since the period around the 6-cycle of the Riemann surface is nonzero, this cycle 
cannot collapse, and the theory cannot have a phase transition. 

For large positive x we have m ~ f and = ^ ~ f • Therefore the two simple zeros 
at z± = iu) ^ 'i\fu ~ ± \fx are approximately on the real axis and are far separated. In 
this limit it is natural to draw the cuts from these points to infinity. 

As X becomes smaller and negative the two zeros move more into the complex plane. 
They move to the upper (lower) half plane for q positive (negative). Let q be positive. As 
X —00 they approach iu ± 2^/u ~ Zi^— | ± y^^, i-e. they move up along the positive 
imaginary axis and approach each other. Here it is more natural to connect the simple 
zeros of y'^ by a cut rather than running the cuts to infinity. This amounts to performing a 
modular transformation on the Riemann surface. In this configuration the period around 
this cycle appears to measure the number of D-branes there. 

This gives us a clear picture of why the phase transition at x = is smoothed out and 
how the fiux (period around the 6-cycle) is continuously connected to D-branes (period 
around the a cycle). 
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Finally, we would like to suggest that the contour integral ( |3.20|) of the force around 
the b cycle can be interpreted as an imaginary energy difference between the two Fermi 
surfaces. In other words, I = iq is the difference in fermi energies between the two sides. 
This puts the parameters q and fx on somewhat similar footing and might suggest a deeper 
interpretation of these phenomena. 

3.2. The relation between FZZT and ZZ branes 

In this subsection, we summarize and slightly extend the discussion in As 



we have argued above, the operator ( |3.1| ) should be interpreted as the insertion of an 
FZZT brane. It was observed in the CFT analysis that one can analytically continue 
the formulae for FZZT branes. In other words, the branes are labelled by a parameter 
s, and two different values of s could give the same boundary cosmological constant fiB- 
We will interpret this as the two sheets of the Riemann surface, in other words we can 
define the operators W {Pj^{z)) and W{P-{z)). It is of interest to find whether there are 
points zq away from the cuts where Fe//(;zo) = 0. A probe eigenvalue located at these 
points will remain stationary, even though it could be unstable. At these points y — and 
R{P^{zo)) = R{P-{zo)). In fact we suggest that the operator 

Z = W{P+{zo)) - WiP-izo)) (3.21) 

creates a ZZ brane at zq. In the classical limit of string theory, the ZZ brane is infinitely 
heavy and it makes sense to think of it as an operator, a deformation of the theory. 
Similarly, the above analysis in terms of Riemann surfaces is also valid in the classical 
limit. Note that since y is zero at zq we can think of zq as the position of an infinitesimally 
short branch cut. Then ( p.21| ) corresponds to integrating R{z) through the cycle that goes 
through this infinitesimal cut and the main cut. In fact, the formula ( |3.21| ) was inspired 
by the formulae in [|| which express the boundary state of the ZZ brane as a difference 
between boundary states of FZZT branes. 

Let us consider, for example, the bosonic string. In that case the FZZT brane param- 
eter s is given by 

cosh %bs = fx b/ ^/p' (3.22) 

where fl is, up to an unimportant factor, the bulk cosmological constant, and /Us is the 
boundary cosmological constant. Then the ZZ brane boundary state can be written aJl 

\D)zz = \Dis+)) FZZT -\D{s+)) FZZT , s±=i{l±b) (3.23) 



We concentrate on the (1, 1) ZZ brane. 
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Note that = 

One can compute the /i and fXB dependence of the FZZT branes for the bosonic string 
with c = 0. We do this by computing the disk one point function for the insertion of a 
bulk cosmological constant, which equals the derivative of the disk with respect to the bulk 



cosmological constant: d^ZpzzT = {Vb)- We compute this using the formulas in ||37|,p8 
Using ( |3.22|) we express the answer in terms of the boundary cosmological constant to find 
that 



df^^D {2iJB - I^B + (3.24) 

which is the expected form of the singular part of the resolvent for the matrix model 
corresponding to bosonic c = 0. This also gives the force on an eigenvalue. We see that 
the force vanishes if = One can check that this is indeed the value of that 

appears for the FZZT in (|3.23| ). Furthermore, as s varies continuously between and s_, 
the value of ^b goes from fx^ through the cut that starts at fiB = — and back to //^. 

Now let us consider the superstring theory. There are two types of FZZT branes 
distinguished by the sign in the supercharge boundary condition. We associate the operator 
( ^■1| ) to the FZZT brane of OB theory with r] = —1 (in the notations of [|T^) with hb = iz. 
For c = and /U > 0, one finds that the disk expectation value for the rj = —1 brane obeys 

13 



d~^D. ^ ^fi + ^il . (3.25) 

This implies 



d^,D.^ ^iB^fi + l^l (3.26) 

which is the expected answer; i.e. it is the same as the resolvent of the two cut model. On 
the other hand for the f] = +1 FZZT brane we find 

D+ ~ ^iBfi (3.27) 

which is analytic in hb- 

We expect that at negative the expectation value of the ry = — 1 brane will be given 
by (|3.27| ). This is as expected of the resolvent of the two cut model for negative fi. 
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4. The m = 2 case — pure supergravity 

In this section we consider in detail the simplest even critical point of the unitary 
matrix model, which corresponds to m = 2 or /c = 1. As discussed above, we remove Ri 
away by a shift of uj, and find the following equations: 

R2 + ^xRo = , rQ2 = . (4.1) 

The equation r02 = ~-f^2 = is a total derivative, which is integrated to H2 = —q. 
Throughout most of the paper we choose the integration constant q to be real. Since 
H2 = —r'^u), q has a simple interpretation as the boost eigenvalue (the "rapidity"). Solving 
for uj and inserting it into ( [1.1|) , we find 

1 1 0^ 

r" - -r^ + -xr + \ = . (4.2) 

2 2 r'^ 

It follows from the Lagrangian 

+-r - -xr + — . (4.3) 



In [30| the two-cut hermitian matrix model was studied using different conventions. 



In (|2.9|) g was taken to be imaginary; therefore, (3 = i9 where ^ is a conventionally defined 



angle with periodicity 27r, see (|2.9|). Therefore 



r"-- + ^-rier = 0. (4.4) 
This equation of motion follows from the Lagrangian 

2^2^84' ^ ' ' 

In Euclidean space the standard way to impose the constraint on the angular momentum 
is to add iqO' to the Lagrangian (the factor i is due to the presence of a single derivative). 
This method reproduces ( ^?2|) . 

Instead, in ||3^ the "angular momentum" was taken to be real. This corresponds 
to imaginary q in (|4.2|): q = il/2, so that 

r" r^ + -xr r = . (4.6) 

2 2 4r3 ^ ^ 
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The real parameter I corresponds to shifting the left "Fermi level" relative to the right 
one, i.e. to having different number of eigenvalues in the two wells. This can be seen by 
computing the eigenvalue distribution from the leading order solution (neglecting the /" 
term). Note that the equation ( [4.6|) comes from the (Euclidean) Lagrangian ( [4.3| ). We see 
that with q = il/2 this action is not bounded below. This implies that we cannot find well 
behaved non-perturbative solutions. More precisely, if we look at solutions of V {r) = 0, 
we see that, for very large x, r ^ ^/x. As we decrease x we find that at a critical value 
of X the solution becomes complex. We can focus on this critical region by taking / oo 
limit and defining the scaled variables 



/l/3 + /-l/15^^ ^^^^2/3^^-2/15^ 



(4.7) 



Inserting this into ( [1.6| ), taking the I oo limit, and rescaling y and u u hy numerical 
factors, we find the Painleve I equation, which is well-known to describe the double-scaling 



limit of a one-cut Hermitian matrix model 42-441: 



^u" -u^ + y = , 



and ajF = |. 



(4.^ 



Thus, the large / limit corresponds to removing the eigenvalues from one side of the po- 
tential, and filling the other side near the top where we recover the single cut critical 
behavior. 

On the other hand, for real q the equation ( [4.2| ) has a smooth solution. We analyze 
these solutions below. 



4-1. Solving the "q- deformed" equation 

Let us examine the solutions of (|4.2| ) in more detail. For large positive x we find 



1 



-73 + 36q^ ^ 10657 - 7048 q^ + 1040 q'^ 



4x5/2 32x11/2 128x17/2 
13912277 + 10303996 q'^ - 2156528 + 144704 q^ 



u{x) = r2(x)/4 



X 



+ q 



1 



2048x23/2 



1 



2x^ 



+ q 



9 



5 (-29 + 4^2) (-83 + 12g2) ^ ^ / 1 



2x 



11 



X 



14 



+ 0{x 



-29/2^ 



2 ^ 14 (g2 - f ) 

5 ^^^^ 



(4.9) 
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and for large negative a;, 



r{x) 



2i/V 






X 


1/4 



^ \q\V2 ^ 5 + 18q2 1^1(107 + 110^2)^ 



4|x|3/2 ' 32|a;|3 
2285 + 13572^2 + 6188^4 



20481x1 



128|x|9/2 
+ 0(|xr^^/2) 



u{x) = r2(x)/4 



\q\V2 g2 
4|x|V2 4|a;|2 



+ 



5|q|(1 + 4q2)^ g2 (7 + 8q2) 



(4.10) 



64|xr/2 



X 



+ 



11 


q 


(105 + 664^2 + 336^4) ^ 


2048 


X 


13/2 



+ 0(|x|-8) . 



The fact that there are terms in the free energy non-analytic in g at g = suggests that 
in the dual string theory there is no R-R vertex operator that corresponds to turning on 
q continuously: had there been a standard vertex operator which couples to q, its n-point 
functions at q = would have been non-singular. This suggests that q is quantized. 

It is possible to argue that the asymptotic expansion ( [4.10| ) matches onto the expansion 
( [1.9|) as follows. The differential equation ([4.2| ) comes from the action 



S 



dx 



V2 + i(r2 
2 8^ 



2 



dx 



V2 + y(r2: 



(4.11) 



This action is bounded below. We can find an approximate variational solution by ne- 
glecting the derivative term and minimizing the potential for each x independently. This 
gives a continuous function of x. For q > the function is smooth; for g = it has a 
discontinuous first derivative at x = 0. Including the second derivative term will lead to 
some changes, especially near x ~ 0, but a solution will exist since it is clear that the 
action has a minimum.0 

Since we started with a well-defined and convergent integral (|2.1|), we should end up 
with a finite and real answer for the free energy F. Therefore, it is natural to expect that 
the differential equation has a unique real and smooth solution. Indeed, the argument 
above shows this. It is important that this is the case both for zero and for nonzero q. 
Note that, in order to select the appropriate solution, it is important to set boundary 
conditions both at x ±oo. 



It is amusing to also present two exact solutions of our equations which do not satisfy our 
boundary conditions. For g = ±i the problem is solved with r = ^/x, n = x/4 and for g = ±| it 
is solved with r = a/^; + ^, n = f + ^. 
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In the worldsheet interpretation of these theories we associate the term of order 
i-3(/i+&/2) -j^ ^ with worldsheets having h handles and h boundaries.! The q depen- 

b 



dence in the expressions (|4.10| ) (|4.9| ) is consistent with this scahng and with a factor of q 



arising from a surface with h boundaries. This explains why the coefficient of x^~^^/'^ is a 
polynomial in q of degree n. 

Note also that in terms of the original rank of the matrix, A^, before taking the double 
scaling limit for positive x we have an expansion with only even powers of while for 
negative x we have an expansion with both odd and even powers. 

Ji..2. A large q limit 

An interesting limit to consider is g — > cxo, x — »• ±00 with t = q~'^/^x fixed. This is the 
't Hooft limit with being the 't Hooft coupling. After we define s = q~^/^r, equation 
( ^72|) becomes 

^s^dfs-s^ + ts^ + 2 = Q. (4.12) 

In the large q limit the first term is negligible, and we end up with a cubic equation for 

v{t) = 

v^-tv^=2. (4.13) 

The solution of this equation leads to a free energy F = q^fit) where v{t)/A = d^fit). 
This is exactly the expected behavior in the large q limit when thought of as a large 
limit. 



For generic t (|4.13|) has three solutions. Only one of them is real for all t (this is easy 
to see for t ^ 0) 



t+(t^ + 27 - 3^81 + 6t3^ ' + + 27 + ^^^l + Qt^^ 



(4.14) 



(here the branches of the two cube roots should be handled with care). For t > —-^r/s 
the arguments of the square roots are positive and v{t) is real. It is clear from the form 
of ([4.14| ) that the half integer powers of t cancel when we do the expansion for t — 00. 
For t < —-^T/j the arguments of the square roots are negative and the second and third 
terms in (|4.14| ) are complex, but v{t) is real. As we move in t and the argument of the 



cubic root moves in the complex plane and it is important to keep track of the branch 



Intepretation of this matrix model in terms of world sheets with boundaries was proposed 



already in p3^ . 
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of the cubic root. When we expand for very large t — > — oo we are are then saying that 
(t^ + ■ ■ ■)^/'^ = e^^^'^/'^t + • ■ ■ for the second or third term in ( f4.14|) . This imphes that there 
is a term where the two square roots add, which leads to the half integer powers of t. The 
expansions of v{t) for large negative and large positive t 

W + ^ W + ^^ + '^^t^)^ (4.15) 

t + ^-l + f -f? + o(^), t>0, 

reproduce the highest powers of q in the expansions ( [4.1UD ( |ir^ after remembering that 
u = v/4. It is interesting that one smooth function (^4.14| ) captures the limiting behavior 
of large q. 

Equation ( [4.13|) is obtained from ( [4.12|) by neglecting the derivative term. In the 
picture of a particle with coordinate s (or r) moving in time t (or x) we presented above, 
this is a limit where we neglect the acceleration term and keep only the potential term. 
The particle is forced to stay at the stationary points of the potential. This suggests that 



we can take (|4.14|) , or more precisely s{t) = ^/f (t), as the starting point of a systematic 
expansion of the solution of ( [4.12|) in powers of l/g^, even though usually such an expansion 
with the highest derivative term is singular. 

We interpret this limit as an 't Hooft limit where only spherical topologies, perhaps 
with boundaries, survive. For q = the theory exhibits the Gross- Witten transition 
at a; = 0. Namely F" = x/4 for x > and F" = for x < 0. It is known that this 
transition can be smoothed by the genus expansion with q = 0. Now we see that it can 
also be smoothed by the expansion in the 't Hooft parameter for infinite q. 

An interpretation of this result is the following. For negative t we have D-branes and 
the power of is the number of boundaries in the spherical worldsheet. For positive t 

the D-branes are replaced by flux; therefore, we have spherical worldsheets with insertion 
of RR fields. Each RR field comes with a power of but their number must be even, 
and hence the expansion in powers of t~^. This is the RR field discussed in section 7. 
Note that the power of t agrees with the KPZ scaling of this operator. It is extremely 
interesting that the theory exhibits, both for finite q and infinite q, a smooth transition 
between the two domains of positive and negative x. This is reminiscent of geometric 
transitions P5| - ^ . The fact that the theory is smooth at t = 0, and the strong coupling 
singularity is smoothed by the RR flux, may have implications for QCD-like theories which 
arise from the conifold with D-branes |]46| - |49[] . 
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5. The m = 4 theory 

Here we briefly discuss the m = 4 theory which has some new phenomena compared 
to the m = 2 case. In particular, for a; < we find that the free energy is discontinuous 
with respect to turning on the parameter q. This happens due to the fact that for g = 
and negative x there exist three different solutions: the trivial one where r = 0, and two 
non-trivial Z2 symmetry breaking solutions with the sphere free energy scaling as |a;|^/^ 
(they are related by the Z2 transformation uj —uj). For q = the trivial symmetric 
solution matches to the positive x solution. For non-vanishing q the positive x solution 
matches one of the non-trivial symmetry breaking solutions. It is clear that the structure 
of solutions gets even richer with increasing m. 

The basic equations for m = 4 are 

i?4-^^ = 0, H, + q = 0. (5.1) 

where i24, are given in appendix C. For x > and q = there is a solution with uj = 0, 
and the free energy is given in ( ^.4| ). Deformation of this solution by q is straightforward, 
and we find 



3 V 3 16/ 3 V 128 9 81 



2q 2q f 80q^ 
uj[x) = h — 



3x 3 V 27 



<fF _ r\x) _ ^ 64g^-15 ^_, ^ /_ 1757 ^ 24V _ 560^ j ^ 



dx^ 4 4 144 V 2304 54 243 

(5.2) 

The free energy contains even powers of q only, so we identify turning on q in the unitary 
matrix model with turning on R-R fiux ~ g in the dual string theory. 

Now let us consider a; < 0. If q = then there is an obvious trivial solution where 
r(x) = 0. We have not been able to find a real deformation of this solution produced by q. 
For such a deformation one expects uj to behave as while r ~ However, 

( ^.1|) require that < 0; therefore, the solution is complex. 

There are two other, less obvious, solutions where both r(x) and uj{x) are non- 
vanishing. One of these solutions has the following asymptotic expansion for x —00: 

r(x) = (2|x|/7) V4 - 2^/4 . 7i - (7/2)^/4^1^ + Odxp^^/^) 

^ ^ ^11/^ 48|x|9/4 ^ ' ' 256|x|i9/4 ^' ' ^ 

4 2v/M 48|x|2 ^''^^ 4608|x|9/2^^' ' ^' 
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and the other is related to it by a; ^ — Since these new solutions break the Z2 symmetry 
spontaneously, in the string interpretation they must involve background RR fields. This 
can be shown explicitly by computing a nonzero one point function of these RR fields using 
the matrix model. The RR fields which condense in the broken symmetry background are 
present in the OB theory but not in the OA theory.! We also trace the fact that the m = 2 
unitary matrix model does not exhibit symmetry breaking to the absence of continuous 
R-R parameters in this model. As m increases, so does the number of R-R fields that can 
condense and produce different solutions. 

Both solutions described above possess real extensions to non-vanishing q. Thus, 
we can expand in small q around the broken symmetry solutions, but not around the 
symmetric solution. This means that it is impossible to vary g in a continuous fashion 
around the symmetric solution. A non-analyticity in q was already obeserved for m = 2, 
but for m = 4 it becomes a more dramatic phenomenon related to the symmetry breaking. 
Just as for m = 2, in the full nonperturbative m = 4 theory there is no vertex operator 
which changes q in the symmetric background. This suggests that q is quantized. 

The requirement of matching an x — cxd expansion to the solution ( |5.2| ) for large 
positive X appears to pick one of the two broken symmetry solutions for q > 0, and the 
other for q < 0. The asymptotic expansion of the q > solution for x —00 is 



r_M-\^ + ^ + - 15 _ ,u 32g^-9 ^.,-18/4. 

4 2yn^x/3 -7^4.(21x1)3/4^ 144|xP '^'/'^ ^96v^|x|i3/4 + "^^"^1 > 

(5.4) 

Remarkably, the torus term in the free energy, In |x|, is the same as for x > (see 

(pT2|)). Perhaps this means that the type OB string theory dual to the m — 4 model has a 
symmetry under a change of sign of x. 

To obtain the solution for g < we act on the above with the transformation q — 
uj —uj. The free energy has the same structure as in the k = 2 complex matrix model 
with positive x that will be presented in the next section, although the coefficients are 
different. The key question is whether the expansion ( |5.4D for large negative x matches 



^ In the next section we will discuss a description of the OA theory by a complex matrix model, 
and we will not find any new solutions exhibiting the symmetry breaking phenomenon. 
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onto the large positive x expansion in ( |5.2|) . We believe that this is the case: a clear 
argument in favor of this can be given in the limit of large q. 

To take such a limit, we assume q > Q and define uj = q^^^h, x = q'^^^t, = q^/^v. 
Then the derivative terms in are suppressed, and we find 

^ o ^ (5.5) 

vh{ — v + 2h'^) = 1 . 

From the second equation we solve for v 



Prom the expansions ( |5.2| ), (|5.4| ) we see that h is negative. For h < we can only choose 
the solution v+. 

Now we substitute this into the first equation in (|5.5| ), move the square roots to one 
side of the equal sign, and then square the equation to remove the square root. After all 
this we find the equation 

-12 - 864/i^ + 448/i^° - 36ht - 96h^t - 27h'^t'^ = . (5.7) 

In order to analyze it, we define the variable 

y = ht. (5.8) 

Then the equation becomes quadratic for 

-12 - 864/i^ + 448/i^° - 36y - 96h^y - 27y'^ = . (5.9) 

Once we find h as a. function of y we can find t = y/h. The solutions of this equation are 



,5 _ 54 + 6?/ ± 5^v/28 + 3{y + 2)^ 

h± - ^ . (5.10) 

The solution with the plus sign, is always positive and nonzero as a function of y. As 
we remarked above, we are not interested in solutions with positive h. The solution h- is 
non-positive, and it is zero only for 

yo = -2/3 . (5.11) 
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As a function of y, /i_ has a maximum at this point. 

As y ^ ±00 we find that ~ c±|y|^/^, where c± is a numerical constant. This 
imphes that t diverges as y^^^ and that h ~ 1^1^^'^- This agrees with the expansion ( |5.4| ) 
which apphes for negative t. If we are interested in a solution where h decreases as t ^ cxd 
it is clear that we should look at the region near yo = —2/3. This forces us to choose the 
solution h^. As we change t continuously, we will stay on the /i_ branch of (|5.10| ) since 



the two branches never cross. For y close to yo we find that h ~ yo/t, which agrees with 
the expansion ( p.2|) for large positive t. 

It is clear that, as t decreases from infinity, y departs from yo. Since is always 
negative, in order to have negative values of t we need to have positive values of y. We 
see then that the relevant region is y > yo- When y — ^ yo we have t 00; when y = 
we have t = 0, and when y — * +00 we have t — > —00. For large positive y, h^_ — — 9y/56, 
which implies 



This agrees with the limiting form of the solution for t —00, (|5.4|) . Expanding the 
solution of ( |5.1U| ),( pT^ further for large negative t, we find 



For large positive t we find the expansion 

The coefficients in these expansions agree with the leading powers for large q in the ex- 
pansions of uj{x) in (p.4|),( p?^ . This shows that an appropriately chosen solution of ( |^ ) 
indeed interpolates between the expansion (|5.2| ) and (|5.4|) in the limit of large q. This also 
suggests that there is a smooth interpolating solution for finite q. Presumably, this could 
be checked through numerical work. 

6. Complex matrix models 

The complex matrix models are based on an (A^ -f x A complex matrix M with 
partition function Z, free energy F and potential V pl| - |35[] 

P (6.1) 

V{MM^) = ^gj{MM^y . 
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Using the U{N + q) x U{N) symmetry we can bring M to the form Mij = XiSij with 
A, > 0. Then (16.11) becomes 



U dy,y!e-Ty(y^^A{y)\ (6-2) 

■ -1 



1=1 

My) = ll{y^-yJ) ■ 

i>j 

Consider first the large limit with a fixed 7 - not the double scaling limit. The 
second form of the integral in ( |6.2| ) in terms of y is similar to the standard hermitian 
matrix model with two exceptions: the factor of y^ in the measure which can be written 
as a contribution to the potential of the form g^^lny^, and the fact that the domain 
of yi is restricted to be positive. Let us now examine the first form of the integral in 
We replace X]'^'^'^ with |Ai|^+^'3', and extend the range of the integral over each A^ to 
(—00, +00) (this adds an inessential constant to F which is independent of the parameters 
in V). Introducing the eigenvalue density p{X)dX the effective potential is 

Veff = dXp{X) (v{X'') - ^^^^^ In |A|^ - dXdX' p{X)p{X') In jA^ - X'^\ 

dXpiX) (viX') - In |A|) - 1 dXdX'piX)piX') In |A - A'| 

(6.3) 

where we have neglected terms of order we used the fact that V is invariant under 
A — > —A and we defined p = p{X) +p(— A) with p(A < 0) = 0. In appendix E we derive the 
loop equation for this model and we analyze the solution. 

For q <^ N ( |6.3|) is exactly like the standard hermitian matrix model with an even 
potential (up to an inessential over all factor of 2). We conclude that in this case the large 
limit of correlation functions are identical to those of even operators in the standard 
hermitian model. Furthermore, if V{X'^) and are such that the complex matrix model 
has one cut that does not reach A = 0, then the corresponding hermitian model would have 
two cuts. If the cut in the complex matrix model reaches A = then the hermitian matrix 
model has just one cut. The important lesson is that the planar correlation functions are 
the same in the two models. This will be important for us since we will interpret these 
two models as the OA/OB version of certain string backgrounds. 
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For q N ( |6.3| ) is like the standard hermitian matrix model with an additive loga- 
rithmic potential. This is the Penner model. We will not explore it further but will assume 
that q is finite when N ^ oo. 

Consider now the corrections and the double scaling limit. We should distinguish 
three classes of critical behavior 

1. A cut of eigenvalues ends at generic A. This critical behavior is described by the 
standard KdV hierarchy of the single cut hermitian matrix model. 

2. Ends of two different cuts meet at generic A. This critical behavior is described by 
the mKdV hierarchy because it is the same as in the two cut hermitian matrix model 
or as in the unitary matrix model. 

3. A cut ends at A ~ 0. In the leading order in (with fixed finite q) this critical 
behavior is the same as in two cut models. However, because of the measure factor 
|Ai|^"'"^'^, the corrections are different than in those models. 

Since only the third kind of behavior is new, let us focus on it. The string equation is 
a differential equation for 

u = 2dlF{z) (6.4) 



where z is a scaled version of A^. The string equation in this case pT|-|3i 



1^101 



- \'^[u]dl'R'[u] + \{d,n[u]f = q' 



1=1 



(6.5) 



(6.6) 



where Qi[u] are the Gel'fand-Dikii differential polynomials 

Qo[u] = ^ 

Q2[u]=-^{dlu-?,u') 

and 7?.[w] is the string equation of the hermitian matrix model with couplings ti. Usually 
the sum over / in 7?. is finite going up to /c. Then we normalize the highest ti such that the 

n = u^ ^ .... 
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For the simplest case, we provide a derivation of this equation in Appendix D. 
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For q = equation (|6.5| ) is satisfied by a solution of tlie differential equation TZ = 0. 
This is exactly the equation of the standard hermitian one-matrix model in its one cut 
phase. The fact that we find this equation can be easily understood by examining the 
expression for the integral (|6.2|) in terms of y. For q = this is exactly the same integral 
as in the one-matrix model, except that the integrals over the eigenvalues are only over 
positive y. As long as the support of the eigenvalue distribution is away from y = 0, the 
effect of the restricted range of integration is nonperturbative. Indeed, for 2; > 0, it is 
possible to check that to all orders in the 1/z expansion ( |0| ) leads to 7^ = 0. Clearly, 
this is not the case for z negative when the effect of the restricted range of integration is 
important. In fact, the nonperturbative corrections to the expansion in 1/^ for z positive 
are also sensitive to the bounded range of y. Therefore, these nonperturbative effects are 
captured by (|6.5D but not by 7^ = 0. 



6.1. The k = 1 case 



For a general q, the string equation of the first nontrivial case /c = 1 is 

un^[u] - ^7^[w]a27^[w] + hd.niu])^ -q^ = o 

TZ[u] — u — z. 



(6.7) 



Using the substitution of [31 



u{z) = f{zf + z 



{<0.i 



the string equation becomes 



dlf -f-zf + l^=0. 



(6.9) 



After rescaling / = 2 ^/^r and it becomes the equation we found in the two-cut 

matrix model, ( [4.2|) , but with the sign of x reversed |31]. Note also that the equations 
(|6.8|) for the free energy become, upon the rescaling. 



d^F 1 



r^(x) 



(6.10) 



up to a non- universal term ~ x. This is consistent with the normalization of the free 
energy in (|2.6|). 
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We conclude that the complex matrix model with positive z is the same as the two cut 
theory with negative x and vice versa, i.e. z = — x.0 Therefore, asymptotic expansions of 
the solution ( [4.9| ),( ^?TI1| ) apply to the complex matrix model as well. The complex matrix 
model with parameter q describes, for positive x, type OA string theory in the background 
of q DO-branes. This agrees with the fact that, in ([4.101 ), the maximum power of q in each 
term in the free energy corresponds to the number of boundaries in a spherical worldsheet. 
Our interpretation of the complex matrix model at negative x, described by the expansion 
{WM) , is in terms of closed string in presence of R-R flux. 



6.2. The k = 2 case 

In the k = 2 case some new phenomena appear compared with k = 1. First, let us 
consider q = 0. For z > 0, the k = 1 complex matrix model is topological, but the k = 2 
complex matrix model has a non-trivial genus expansion. Indeed, for k = 2, the equation 
is H 

w7^2 - 7^7^'72 + (7^')V4 = , n = u'^ -u"/3- z . (e.ii) 

Clearly, it is solved by any solution of 7?. = 0, which for /c = 2 is the Painleve I equation. 
Around z ^ oo the solution of this equation has the asymptotic expansion 

1 49 

but it cannot be continued in a smooth way to negative z. 

Luckily, there are solutions of (|6.11|) which do not satisfy TZ = [|5^. Neglecting the 
derivatives of u we see that we can solve ( |6.11| ) for negative z without satisfying 7^ = by 
u = 0. Expanding around it we find 

1 225 6906075 



''<^' = -4iJ-32l7- 1024lI5+O(.-"). (6.13) 



The authors of [|5^ showed numerically that the two asymptotic expansions ( |6.12| ) and 
( |6.13| ) are connected through a real and smooth solution of ( |6.11| ). 

Generalizing ( |6.11|) to include q, we find (|6.5|). Now the asymptotic expansions of the 



solution are, for positive z, 



u{z) = + J4 - + 0{z-'''% (6.14) 

2zi 24 2:^ 



The boundary conditions are u z^l/m) (0) for z — > +oo (— oo). We also demand that TZ 
is positive and goes to zero as z ^ +oo. See appendix D for further explanation. 



27 



and for negative x, 

u{z) = {cP - ^) + " " + (^(^-")) . (6.15) 

The g ^ limit of this solution is ( |6.13| ). Thus, in the k = 2 complex matrix model we do 
not find a discontinuity at g = that is present in the m = 4 unitary matrix model even 
though the two models have the same critical exponents. The discontinuity of the unitary 
matrix model is related to the Z2 symmetry breaking. From the string point of view, the 
RR fields that condense in the OB string are absent from the spectrum of the OA string 
dual to the complex matrix modeL^ 

We identify this theory with the OA version of the (2, 8) superminimal model for 
which the KPZ scaling of the h handle and b boundary surface is 2;~^(''+^/^~^)/^~^. The 
expressions above are consistent with this identification, with the power of q being the 
number of boundaries or insertions of the RR ground state operator. 

As for /c = 1, on one side we have boundaries - D-branes; on the other side we have 
no boundaries but only fiuxes, hence only even powers of q appear. 



7. String theory in one dimension 

Here we consider the type theory in one target space dimension. This theory is also 
known as pure supergravity or c = noncritical string theory. Using standard Liouville 
conventions we have Q = ^ind 7 = Therefore, the contribution of worldsheets with 
h handles and b boundaries scale like ^3(2-2/i-b)/2 

The theory has a single closed-string NS operator, the worldsheet cosmological con- 
stant. In the (—1, —1) picture it is e^'^. For ^ = 0, using free field theory, we find two 
candidate R-R operators. In the (—1/2,-1/2) picture they are given by V± = a^e^'^l'^, 
where cr^ are the two spin fields of the super-Liouville theory. The worldsheet (— 1)"^ 
projection in the OA theory leaves only V"_j_, and in OB it leaves only V-. For nonzero /i 
their wave functions are determined by solving the minisuperspace Schrodinger equations, 
as in [1^]. Because of the behavior at (/> ^ cxd, only one of the two operators is acceptable. 



Our conventions are such that for < only is acceptable, and for (U > only V- is 



It is amusing to point out three simple, but unphysical solutions of ( [6.51 ). For g = 1/2 it is 
solved by u = 0, for g = 3/2 it is solved by « = 2/2;^ and for g = 5/2byn = 6/2;'^. More generally, 
it is easy to show that the solution can have double poles with residues 2 or 6. 
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acceptable. This means that in the OA theory there is one R-R operator V+ for n < and 
no R-R operator for fj, > 0. In the OB theory the situation is reversed: there is V- for 
H > and no R-R operator for fx < 0. 

More generaUy, reversing the sign of n can be undone by transforming the worldsheet 
LiouviUe fermions as i/^ — — V' with if) ^ Since these are the only fermions in the 
problem, this has the effect of changing the sign of the worldsheet (—1)^ in the R-R 
sector .111 Therefore, it amounts to reversing the projection in this sector and interchanging 
OA and OB. We conclude that the OA theory with fi is the same as the OB theory with —fi. 
Unlike the situation in c = 1 where both theories are invariant under fi —fi [|T7| , here 
they are interchanged. 

All this can be summarized by a simple "spacetime" picture.0 Recalling the situation 
in the c = 1 theory, and reducing it to one dimension, we find only one field. In the OB 
theory it is C; in the OA theory it a component of a gauge field At in what was originally 
the Euclidean time direction. Since OA and OB are related hj n ^ —fx we can focus, 
without loss of generality on ObIiI. The leading order term in the action for C is 

<ple'^(^\d,Cf (7.1) 
where T{(f)) = fxe^'^ is the tachyon field. The equation of motion is solved by 

C{<p)= re-2^(^')#'. (7.2) 



For fx > this solution is well behaved as ^ ^ oo, and leads to a "fiuctuating degree of 
freedom." Since it is linearly divergent as ^ ^ — oo, it is not a normalizable mode and 
corresponds to a physical vertex operator |]5I||T7| ] V-. For /i < this solution is badly 
behaved for (/> ^ cxd and therefore should be discarded. 

We can introduce background charge q at infinity which leads to background field 

C{<l)) = iq /'^e-2^('^')#'. (7.3) 



This fact is familiar from the study of the Ising model. Changing the sign of the fermion 
mass (moving from one side of the transition to the other), the order and the disorder operators 
change roles. This means that we change the sign of the ( — 1)^ projection in the R-R sector. 

We put the word spacetime in quotation marks because we have only space parametrized by 

4>. 

There is also a RR one form A^,. 
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This can be done by adding the surface term 

iq J d(t)'d^:C{(f)') = iq (C(oo) - C(-oo)) (7.4) 



to (|7.1|) . (We view the action ( [7.1|) as a Euchdean action and therefore we put an i in 



front of the topological term d(j,C. This leads to imaginary classical solution for C{(p) as 
is common in classical solutions of D-instantons.) For // > 0, integrating out C{(p) we find 

2 

^ log(A/|/x|) where A is a cutoff on cj). 

For fx < 0, the solution of ( |77^ ) is badly behaved. In this case it seems possible to 
add charged ZZ branes that will absorb the flux of C and somehow lead to a finite answer. 
These are charged D{—1) branes and their number should be precisely equal to the flux. 
Therefore, from the string theory point of view it seems that q should be quantized. 

More generally, it is natural to assume that q arises from the charges which exist in 
the theory. These are ZZ branes, which are localized near (p = oo. Therefore, q should be 
quantized. With integer q the surface term ( [7.4|) does not change if C(oo) C(oo) + 27r, 
as expected from the periodicity of C. 

q appears as a ^ term, but unlike similar terms, the theory is not periodic in q because 
the standard process of pair creation which can screen its charge does not take place 
because the D-branes are locked at infinity. These D-branes look like instantons. However, 
unlike with ordinary instantons, we do not sum over q but keep it fixed. There are two 
reasons for that. First, the standard argument which forces us to sum over q relies on 
creating an instanton anti-instanton pair and separating them. This argument does not 
apply here because the D-branes are forced to be at one point in our "spacetime." Second, 
backgrounds with different values of q differ by infinite action. 



8. R-R physical vertex operators 

In this section we classify the R-R vertex operators, which distinguish the OA theory 
from the OB theory. 

Consider a "matter" SCFT with central charge c and an R-R primary operator with 
dimension A > Before performing any GSO projection there are two such operators 
corresponding to the states |±)m, where the sign denotes (— 1)"^ with F the worldsheet 
fermion number. It is convenient to consider the sum and the difference of the zero modes 
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of the left and right moving matter supercharges, Gm = -^{GM,Left + ^Gm, Right) and 
Gm = -^{GM,Left - ^Gm, Right)- They act of these two states as 



Gm 


+)m 


= v 




^ 1 \ 


Gm 


—)m 


= 






Gm 


—)m 


= V 


/a- 




Gm 


+)m 


= 0. 







There can also be other representations in which Gm\+)m = Gm\—)m = 0, but Gm\—)m 
and Gm\—)m do not vanish. Such representations are not present in the superminimal 
models. But they do exist in more generic systems. For example, they exist in c = 1 and in 
the flat space critical theory, where we can think of the "matter part" as the superconformal 
field theory of nine free superfields. We will focus on models where all the representations 
are of the type (|8.1|). 

We couple this system to Liouville with central charge cl = 10 — c. Physical vertex 
operators in the R-R sector have ghosts. In the {—\:~\) picture the total dimension of the 
matter and Liouville is |. Therefore, the dimension of the Liouville operator needed 
to dress the matter operators of dimension A satisfies A^, — f|- = — (A — We denote 
the two R-R operators with these dimensions as |±)l- (Since A^ — ^ < 0, these are 



operators and not normalizable states fSTI]; nevertheless, we will use the state notation.) 
The action of the zero modes of the Liouville supercharges Gl and Gl (which are again 
linear combinations of left and right moving supercharges) on them are 



Gl\+)l 


= V 


Gl\—)l 


= 


Gl\-)l 


= V 


Gl\+)l 


= 0. 



(8.2) 



The total superchargest^l are G = Gm + Gl and G = Ga/ + Gl of (|T|) and (|3). 
In the OA theory the candidate operators in the { — \^~'k) picture have (—1)^ = 1, and 



Actually, in order to obey the proper anticommutation relations we need a cocycle. The 
proper expression is G = Gm x 1_l + (—1)^*^ x Gl where (—1)^*^ is the matter fermion number. 
We have a similar expression for G. In order not to clutter the equations, we will suppress the 
cocycle. 
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therefore the aUowed operators are \+)m\+)l and \—)m\—)l- Imposing that they are 
annihilated by the total supercharges we find that there are no such physical states. 

In the OB theory the candidate operators in the (— ^5 — |) picture have (—1)^ = — 1, 
and therefore the allowed operators are \+)m\—)l and |— )m|+)l- Imposing that they are 



annihilated by G = Gm + Gl and G = Gm + Gl of (|8.1|) and ( |8.2|) we find one physical 
operator 

\+)M\-)L + iHM\+)L. (8.3) 

As we mentioned above, in theories with c > 1, there can be matter operators which 
satisfy ( |8.1|) with the interchange of + ^ — . For such operators, the above conclusions 
about the spectra in OA and OB are interchanged. 

Consider now the R-R ground states which exists in those SCFT which do not break 
supersymmetry. It has A = Here there is only one state \+)m satisfying 

G\+)m = G\+)m = 0, (8.4) 

which we take to have even fermion number by convention. In some examples like the 
c = 1 theory there are two such ground states | ±) m and the discussion below is modified 
appropriately. 

In the free field description of Liouville there are two states with Al = ^ which 
satisfy 

Gl\±)l = Gl\±)l = 0. (8.5) 



For nonzero only one of them is well behaved as (/> ^ 00 |T^. For < it is 1+)^ and 



for > it is Imposing the GSO projection we find the physical operators. The OA 

theory has an R-R ground state operator (in the (— I, — ^) picture) 

\+)m\+)l for n<0 in OA (8.6) 

and no such operator for // > 0, while the OB theory has no R-R ground state operator for 
/X < 0, and 

\+)m\-)l for n>0 in OB. (8.7) 

The spectrum of NS operators in these theories is independent of the sign of fx, and it 
is the same in the OA and OB theories. Let us summarize the spectrum of R-R operators. 
In the OA theory there are no R-R operators for fi > and for fj, < only the R-R ground 
state is present. In the OB theory with fj, > there is a physical vertex operator for each 
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R-R primary. For n < the spectrum is the same except that the R-R ground state is 
absent. 

One interesting special case is c = 0. Here the general picture above is valid, but 
there are no R-R operators other than the ground state. Therefore, the spectrum of the 
OA theory with one sign of n is the same as the spectrum of the OB theory with the opposite 
sign of This is consistent with our general claim that fi — > —fi exchanges OA and OB in 
this case. 

Another interesting case is c = 1. We can also apply the above construction, keeping 
in mind that there will be two types of matter operators, the ones in and similar 

ones with plus and minus interchanged. There are two R-R ground states in the matter 
theory with opposite fermion numbers. For every sign of n one of them leads to a physical 
vertex operator in the theory coupled to gravity. This is consistent with the fact that these 



theories are invariant under n [[17 



In nonunitary theories the lowest dimension matter operator is typically not the 
identity operator. Then, a generic perturbation of the SCFT is given by this operator 
X f d^OOmin, and upon coupling to Liouville, this operator is dressed. We assume that in 
the conclusions ( p^) (|8.7| ) about the existence of the R-R ground state, we simply have to 
replace fx by the coefficient x in this case. 

The R-R ground state operators ( |8.6|) and ( p77|) represent fluxes. If these fluxes can 
only be induced by D-branes, they are quantized and cannot be changed in a continuous 
fashion. Then, these operators can appear as vertex operators in the perturbative string 
theory, and they can appear with quantized coefficients in the worldsheet theory. But they 
cannot exist as standard vertex operators in the complete nonperturbative theory. 



9. Superminimal Models Coupled to Supergravity 

In this section we review some basic facts about (Ap_i, Ap'-i) superconformal minimal 
models and their coupling to super-Liouville theory. 

The central charges of the superminimal models are 

c = 1 . 

pp' 

The operators are labelled by two positive integers j and j' subject to 

1</<P-1, 1<J<P-1, jp'>j'p- (9.1) 
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Their dimensions are 

8pp' 32 ^ ' 

Operators with even j' — j are NS operators and those with odd j' — j are R operators. 
The R operator with jp' = j'p is the R ground state. It has h = c/16. 

The super-minimal models are characterized by two positive integers p' and p subject 
to restrictions p' > p; p' — p = mod 2; if both are odd, they are coprime, and if both are 
even, then p/2 and p' /2 are coprime. There is also a standard restriction that if p and p' 
are even then {p — p')/2 is odd (see Appendix B). For models with odd p,p' there is no 
R state with h = c/16; hence these models break supersymmetry. For models with even 
p' there is such a state, but we will show that upon coupling to super-Liouville theory, 
it does not give rise to a local operator. 

When coupled to super-Liouville we need the total c = 10. This fixes 



Q=Jl^ = l±4. (9.3) 

V Z y/pp 

The operator labelled by is dressed by e^""^ with exponent 

Zypp 

In the super-Liouville action we may include the coupling to the lowest dimension operator 
in the NS sector of the matter theory, Omin, 

Pmin = „ ^ , P, P even . 

lypp 

In models with odd p, p' the lowest dimension operator is in the R sector, and then we find 

n _ p+p'-l 

A^min „ / — 7 •' — ' 



Since models with odd {p,p') break supersymmetry, one might suspect that their coupling 
to supergravity leads to theories which are equivalent to some bosonic minimal models coupled 
to gravity. Indeed, there is evidence that they are equivalent to bosonic {p,p') minimal models 
coupled to gravity, with the dressed lowest dimension operator turned on. In these bosonic models, 
Q = ^^P±l'^ and /3:^i„ = £±£=1, so that = ^^1. This is the same as in odd ip,p') 

-I /pp' V 2PP' Pmin P+P J- 



superminimal models coupled to super-Liouville and perturbed by the dressed lowest dimension 
operator from the R sector. Thus, the scaling of the partition function is the same. So are the 
gravitational dimensions of dressed operators. 
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Consider now the series of theories {p = 2,p' = 2m) with m = 1,2, .... The standard 
restriction of ||5^ (see Appendix B) requires that m is even. If we denote m = 2k, then 
the resulting (2, 4k) theories, when coupled to super-Liouville, will match in the OB case 
with the critical points of the two-cut matrix models (these critical points are also labelled 
by positive integer k and belong to the mKdV hierarchy found in [25||). 



Theories with odd m can be obtained by starting with theories with higher even m 
and flowing down by turning on a Z2 odd operator. In terms of the worldsheet description 
this is a superconformal field theory coupled to supergravity which is perturbed by a R-R 
operator. Such theories are not expected to be massive field theories coupled to Liouville. 
The reason is that the R-R vertex operator involves the ghosts and it is no longer true 
that the ghosts and the matter fields are decoupled. This is the origin of the difficulty 
of describing strings in background R-R fields in the RNS formalism. However, in the 
matrix model there is no such difficulty and the ZS hierarchy allows us to describe such 
backgrounds. In particular, we find new critical points, the ones with odd m, by turning 
on such R-R operators. For a discussion of two-cut matrix models with odd m > 1, see 
Appendix C. 

The super-Liouville theory is characterized by 

^ m + 1 , . 

Qm = . 9.6 



The operators are labelled by (j = 1,/ = 1, 2, m — 1). The operators with odd j' are 
from the NS sector, and those with even j' are from the R sector.lii The Liouville dressings 
of these operators are determined by 

= (9.7) 
2\/m 



The lowest dimension operator is the operator with j' = m — 1 and 



Pmin — f^m-l — (9.8) 

Since this is the most relevant operator, we will turn it on as the generic perturbation 
with coefficient x. For even m it is an NS operator, while for odd m it is an R operator. 



IS We removed from the list the operator labelled by j' = m. For the SCFT with even m this 
is the R ground state which was discussed in the previous section. 
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Therefore, for even m the theory has Z2 symmetry under which all the R fields are odd, 
while for odd m this symmetry is broken. 

Standard KPZ/DDK scaling shows that the correlation functions on a surface with h 
handles scale like 



<3j('i-i) + ^. 



0,1 



^n.in =ji-'^-^){h-^)-Y.J-4^, (9.9) 



This matches the discussion of operator dimensions that we discussed after ( p.20|) . The 



enumeration of operators, the scaling in (|9.9| ) and the Z2 assignments agree with that in 
the double cut (unitary) matrix model pB|-PB| with the string action 



dxJ2tiHi+i] tm = l, to^x. (9.10) 

1=0 

Here tm-i is redundant and can be shifted away by an appropriate redefinition of the 
couplings. As in it can be interpreted as the boundary cosmological constant which is 



important only when worldsheet boundaries are present. More precisely, it is the boundary 
cosmological constant in presence of the FZZT branes [|37 , |3H] . ti is the coefficient of the 



operator labelled by j' = m — / — 1 in the CFT. For example, ti=o = x is the coefficient of 
the lowest dimension operator labelled by = m — 1, ti=rn-2 is the cosmological constant 
(the dressed identity operator) labelled by j' = 1. 

It is simple to repeat this for the OA theory. The complex matrix model is a Z2 
orbifold of the two cut model. It has the same spectrum of Z2 invariant operators as its 
parent two cut model and the Z2 odd operators are absent. The same conclusion about 
the spectrum applies to the OA theory which is an orbifold of the OB theory. Therefore, 
the OA theory has the same spectrum as the complex matrix model. 



10. The Torus Path Integral 

In this section we compare the matrix model results for OB and OA theories with 
the genus-1 path integral for supergravity coupled to the {Ap-i, Ap'-i) superconformal 



minimal models. For even spin structures this was calculated in ||53|| . 

As in the bosonic case, we perform the integral over (po first, and it contributes a 
volume factor Vl = —(In \x\)/(3min- Integrating over the rest of the modes, we find 

(r,s) 
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where Z'^^ij^ f ) is the matter partition function in the spin structure (r, s). One can repre- 
sent the partition functions of the superconformal minimal models as linear combinations 
of the partition functions of a compactified scalar superfield [^. For the supersymmetric 



(Ap_i, Ayi-\) series, in each even spin structure there are relations [q2| 



Zp,p. = -z(7^/v^ = ^fwTi) - i,z{n/V^ = v^^Tv) , odd 

{ I (10.2) 

Zp,p. = -Z{R/^ = ^JWm - ^Z{n/^/^ = V^TV) , P.P' even , 

where Z is the partition function of the "circle" c = 1 theory, and Z is the partition 
function of the super-afRne theory. Substituting these expressions into eq. ( |1U.1| ), and 



performing the integrals, we find the sum of the contributions of the (— , — ), (— , +) and 
(+, — ) spin structures, each weighted with the factor l/2:lll 



1 {p-rw-i) 

Zeven = ; ^ — In x , », » odd 

16 {p + p' -1) ' ' ' 

1 {p-l){p' + l 

Zeven = "T^ 7 — \ — ; m X , p,p eveu 

16 [p + p — 2) 



;io.3) 



Note that the answers depend on In Our discussion in section 8 about the spectrum of 
these theories shows that the even spin structures are independent of the sign of x. This 
sign is important in the odd spin structure. 

The continuum calculation of the odd spin structure is more difficult and is beyond the 
scope of this paper. However, even without calculating it, we can carry out an interesting 
consistency check of our OB and OA conjectures. Since the odd spin structure contributes 
with opposite signs in the OA and OB theories, we can check that the average of their 
torus partition sums agrees with (|10.3|) . Furthermore, since the even spin structures are 



independent of the sign of x, we should find 

1 



2 ^ZAix) + Zeix)) = Zeven(|x|) . (10.4) 



In models with odd {p,p'), x is the coefficient of the dressed lowest dimension Ramond 
operator. Note that the sum over even spin structures is the same as for bosonic {p,p') models 



coupled to gravity, up to overall coefficient (in the bosonic case the coefficient is 1/24 ||5^). It 
is possible that inclusion of the odd spin structure will give exact agreement with the bosonic 
theory. This would provide further evidence that the odd {p,p') superminimal models coupled to 
super-Liouville are equivalent to {p,p') minimal models coupled to Liouville. 
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Let us consider theories with p = 2,p' = 4k where the 1-matrix models discussed in 
sections 2-5 and in section 6 describe OB and OA theories respectively. In this case the 
worldsheet computations lead to 

^even = -:^ln|x| , (10.5) 

lb 

independent of k. Let us compare this with the matrix models. We start with the complex 
matrix models discussed in section 5, changing the notation z ^ x. If we set g = in 
( |6?5|) , then the perturbative solution for x > is obtained simply by setting 7l{x) = 0. 
This is exactly the set of KdV equations that describe one-cut Hermitian matrix model, 
except the result is divided by 2 to eliminate the doubling of the free energy that appears 
for symmetric potentials: this is the origin of the factor 2 in ( |5.4| ). It follows that the torus 
path integral in the k-th multicritical complex matrix model for x > is the same as in 



gravity coupled to the (2, 2k — 1) minimal model which was calculated in [|5^. Substituting 
p = 2, p' = 2k — 1 into the result of we find 



^P~^^^P'~^yinx = -^lnx. (10.6) 



24(p + p'-l) 24k 

Of course, this result can be obtained also by using the KdV equations 7l{x) = and 
integrating ( |6.4] ). 

For X < we return to ( |6.5| ) and set q = 0. It is easy to see that the asymptotic 
solution of the equations is u ~ ^ (note the agreement with our /c = 1,2 

expressions). Using ( |6.4| ) we learn that 

ZA = l~w]'^i^ (10.7) 

\ -| In |x| X < . ^ ^ 

In the multicritical two-cut models ( |2.4| ), we find that the torus path integral for OB theories 
is 

zb = <:-^i^^ (10. 



X < 0. 

The vanishing result for x < follows from the fact that, for q = 0, the relevant solution 
for large negative x is the trivial one, with r = up to non-perturbative corrections. 
Using ([T077D ,([T0l 



^even = ^ {Za{x) + Zb{x)) = -^ln|x| (10.9) 
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in agreement with the worldsheet value ( |10.5| ). As expected, it is independent of the sign 
of X. Using the matrix model results we get a prediction for the odd spin structure 

1 I ^4Sk ^^-^ X > 

ZoM = - {ZAix) - Zb{x)) = I (10.10) 



2 



In a; < 



which does depend on the sign of x. It will be interesting to check this with an explicit 
continuum calculation. 



11. Conclusions and Future Directions 

In this paper we studied two types of matrix models: 
1. The unitary matrix model (or, equivalently, the two-cut Hermitian matrix model). In 
the double scaling limit it is described by the differential equations ( |2.20| ) 



dr{x) 
5 

du!{x 



/m m 
dxC^tiHi+i + qu) = -^ti{l + l)Ri = 

1=0 1=0 

/m m 
dx(J2tiHi+i + qu^) = J2^i{l + l)Hi + q ^ 
1 r\ 1 r\ 



1=0 1=0 



We identify this model with type OB string theory in a background characterized by 
ti and q. 

2. The complex matrix model, which is described in the double-scaling limit by the 
differential equation ( |6.5D 

uTZ^iu] - hl[u]dln[u] + ^{d^n[u]f = q^ 

oo ^ (11.2) 
n[u] = Y,{i + -)tiQi[u]-z. 
1=1 

We identify this model with type OA string theory in a background characterized by 
ti and q. 

Unlike in the standard Hermitian matrix model, the potentials of these models are 
bounded from below, and we expand around their global minima. Correspondingly, the 
differential equations ( |1 1 . Ij ) ( PT!^ ) have smooth and real solutions. 

We pointed out that, when a superconformal field theory is coupled to super-Liouville 
theory to make a superstring theory, there are in general four independent weak coupling 
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limits. First, we have a two-fold ambiguity in the sign of the odd spin structures. This leads 
to the type OA and type OB theories. Second, we can change the sign of the cosmological 
constant fj, —fi, or more generally, the sign of the coefficient of the lowest dimension 
operator x —x. In the simplest case of pure supergravity (c = 0), these theories are 
identical pairwise: OA with /i is equivalent to OB with — /x. This is not the case, however, 
for more general theories. 

The recent advances in non-critical string theory have been based on the idea that the 
matrix model is a theory of a large number of ZZ D-branes [^[|^,|16|,|1^ . The identification of 
the FZZT branes with the integral of the matrix model resolvent leads to another insight. 
The boundary cosmological constant on this brane can be analytically continued to take 
values on a Riemann surface which is a double cover of the complex plane. The cuts in the 
complex plane represent the eigenvalues; the discontinuity of the resolvent (the derivative 
of the FZZT brane with respect to the boundary cosmological constant) is their density. 
This works well for one of the FZZT branes of the type theory. It would be nice to 
identify the other FZZT brane in the matrix model. 

One novelty of the matrix models is that they allow us to analyze theories with back- 
ground Ramond-Ramond fields. This is an important topic that we have only started to 
analyze in detail. Our preliminary investigation has already led to the following observa- 
tions: 

1. Even without turning on the background R-R couplings ti, the OB theory has solutions, 
e.g. ( |5.3|) , which break the Z2 symmetry that acts as —1 on all the R-R fields. It would 
be nice to clarify the nature of these solutions in more detail§. 

2. By turning on background R-R fields, the Z2 odd ti couplings, we find new theories. 
They involve nontrivial couplings of fields and ghosts, but their description in the 
matrix model is as easy as the even ti fiows. Among them we find scaling solutions 
which are described by the odd m equations in the ZS hierarchy. Unlike the even 
m scaling solutions, these are not super conformal matter field theories coupled to 
supergravity. 

3. The parameter g is a particular Z2 odd coupling. In the OB matrix model it appears as 
an integration constant. In the OA matrix model it can be introduced by considering 
a model of rectangular matrices. Alternatively, it can be introduced by changing 



On the matrix model side, they seem related to performing the matrix integral along other 
contours. 
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the measure in the complex matrix model of square matrices by adding a factor of 
(detMMt)l9l to the measure. In some cases q appears to be related to adding D- 
branes to the system. This is particularly clear in the complex matrix model with 
rectangular matrices. From this point of view it is natural to assume that q must 
be quantized. However, this conclusion might be misleading. If q is introduced by 
changing the measure of the OA model, or if it is introduced as an integration constant 
in the equations of the OB model, we see no reason why it should be quantized. It 
would be nice to find a clear argument which determines whether q should be quantized 
or not. If q is not quantized, we have seen some physical amplitudes which are not 
analytic in q around g = 0. Therefore q cannot vary in a continuous fashion and its 
change is not described by a standard vertex operator. 

4. It is common that for one sign of x the parameter q appears with boundaries on the 
worldsheets and can be interpreted as associated with D-branes, while for the opposite 
sign of X only even powers of q appear, and it can be interpreted as background R-R 
flux. This difference in the behavior as x changes sign is consistent with the behavior 
of the proflle of the R-R fleld strength as a function of the Liouville field What is 
surprising is that the system smoothly interpolates between the behavior at positive 
X and at negative x. Such a smooth interpolation between D-branes and fiuxes is 
reminiscent of geometric transitions p5| - ^ . Our solvable models provide simple, 
tractable and explicit examples of this phenomenon. 

5. By examining the planar limit with x positive and x negative one often finds a phase 
transition associated with nonanalytic behavior of the free energy F . We have seen 
that the finite x (higher genus) corrections smooth out these transitions. Alternatively, 
we can keep |x| large, i.e. continue to focus only on spherical worldsheet topologies, but 
smooth out the transition by turning on q. In other words, the transition is smoothed 
out either by including worldsheet handles, or by including worldsheet boundaries or 
nontrivial R-R backgrounds. These latter cases provide particularly simple examples 
of the interpolation from D-branes to fiux. They can be seen with spherical worldsheets 
and are described by polynomial equations rather than differential equations. In the 
OB theory we have interpreted g as a certain period around a cycle in a Riemann 
surface. The transition associated with the collapse of this cycle is prevented by a 
nonzero period. 

We found many situations in which amplitudes vanish without a simple worldsheet 
or spacetime explanation. For example, the perturbative expansion of F with negative 
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X and q = vanishes for some of our solutions of the OB theory; see e.g. ( 4.10|) but not 
( ^.4|) . Also, the even genus amplitudes in the expansion ( |6.15|) of the OA theory vanish, 
but the odd genus amplitudes are nonzero. Perhaps these vanishing amplitudes reflect a 
deep structure of these theories. 

Interesting insights into the bosonic string counterpart of these systems has been 



gained by interpreting them as topological theories [q5|. It is possible that a similar 



topological structure underlies our examples and our results. For a possible starting point 
for investigating this question, see |j2^ . 

An obvious generalization of the dualities we studied is provided by multi-matrix 
model versions of them. It is likely that they can provide realizations of all {p,p') super- 
minimal models coupled to supergravity. We expect these systems to exhibit a rich struc- 
ture which generalizes the phenomena seen in this paper and in the analysis of the bosonic 
noncritical string. 

We have related the FZZT brane of the OB theory with rj = —1 with the resolvent 
of the two cut matrix model. The FZZT brane with ry = 1 related to the resolvent of the 
complex matrix model. It would be nice to have a description of the two branes within 
the same theory. It would also be nice to have a clearer description of the relation of the 
ZZ branes to the FZZT branes. 
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Appendix A. A simple comparison — the m — A model 

In this appendix we compute some simple tree level correlation functions of RR vertex 
operators using the matrix model results. 
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The pure supergravity example (m = 2) discussed above is particularly simple in 
that the model has only one coupling constant which corresponds to turning on the super- 
Liouville superpotential e'^ from the NS sector. As we discussed in the previous subsection, 
the R ground state is not present as a standard deformation of the theory. Another would 
be R-operator is redundant and can be shifted away. 

For this reason let us consider the next even critical point, m = 4 which we expect 
to correspond to the (2, 8) super-minimal model coupled to gravity. In this model the 
effective Lagrangian is given by 

L = xHi+ eH2 + £Hs - . (A.l) 

15 

Now there are 2 NS operators, corresponding to coupling constants x and x, and one RR 
operator corresponding to e. The coupling x corresponds to the dressed j = 1,/ = 3 
operator of negative dimension, while the coupling x to the dressed identity. The coupling 
e corresponds to dressed R-R operator Or with j = = 2. We will set 5; = and 
calculate the expansion of the free energy in powers of e. 

To perform the sphere calculation we neglect derivatives of r and uj. From the KPZ 
scaling we expect that {OrOr) ~ xlnx, but we do not find such a term. It follows that 
the two-point function vanishes, or is given by a non- universal term ~ x. However, the 
four point function of Or scales as in agreement with the KPZ scaling. 



Appendix B. Super-minimal models with even p and p' 

The super-minimal models are characterized by two positive integers p' and p subject 
to: p' > p, p' — p = mod 2, if both are odd, they are coprime, and if both are even, 
then p/2 and p' /2 are coprime. There is often also a restriction that if p and p' are even 
then {p — p')/2 is odd. The purpose of this appendix is to review parts of the discussion of 
|52} j emphasizing why this requirement is needed. Along the way we will also review some 
useful facts about the super-minimal models. 

The central charge of the super-minimal model labelled by (p, p') is 

c=l-2^P^. (B,l) 

PP 

The operators are labelled by two positive integers m and m' subject to 

1 < m' < — 1, 1 < m < p — 1, mp' > m'p. (B.2) 
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Their dimensions are 

{mp' — va'pY — {p — p'Y 1 — (— 
/iw = + . (B.3) 

Operators with even m! — m are NS operators and those with odd m! — m are R operators. 
The R operator with mp' = m'p is the R ground state. It has h = c/16. 

For even p and p' and odd {p — p')/2 the R ground state is the operator labeUed by 
(m = |,m' = y). For even p and p' and even (p — p')/2 the operator (m = |, m' = y) 
is in the NS sector and it has h = — ^ / = This means that the effective central 

8pp' 16 



charge of the theory |51,56|, which is given in terms of the lowest dimension operator is 
Ceff = c — IGhmin = 1, and therefore the density of states of such a theory is as in the 
c = 1 theory. We conclude that the theory must have an infinite number of super- Virasoro 
primaries, and it cannot be a super-minimal model. Since this argument depends on the 
modular invariance of the partition function, let us examine it in more detail. 



The superconformal characters in the different sectors are [52 



X^ir) = lY.Q^o-^ = xfii,.=o(r) [Kxir) - K~,{r)] 

xrir) = Tt,q^°-^e(^-lf = x^ir + 1) = xg,,=o(r) [Kxir + 1) - K-,(t + 1)] 



Xx 



t) = Trxq"-'-^^ = Xf=i 1 (r) [Kx{r) - K-^{r)\ . 



(B.4) 

Here the traces are in the representation labelled by A = mp' — m'p and we use 

X^ depending on whether m' — m is even (NS) or odd (R). In ( p.4| ) we use the notation 

A = mp' -\- m'p. Xc=i,h are the characters in c = 1 in the same spin structure and 

OO , 9 

(2pp'n+A)^ 

Kx= ^ (^•^) 

n= — oo 

(we moved a factor of r7(r) from K to the factor that multiplies it relative to [|52|). 

As a simple consistency checks examine the leading behavior of the different characters 
as q = e^'^*'^ —>■ 0. Since by ( |B.2D < A < pp' , the leading term in the sum in ( [B.5[ ) is with 
n = 0, and therefore 



2 



Xx [^) '^'^ = Q 



Xxi'^) ^ (l^"" ^''"^spf' = q^~ 



16 
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Under r r+1 each term in the sum ( |.B.5| ) is multipUed by e^^^^^ ^ +("ip p)"-+^^)_ 
If p and p' are even, all the terms have the same phase and K transforms by a phase. If 
p and p' are odd, and m — m' is odd (R representations) again all the terms have the 
same phase and K transforms by a phase. Finally for p and p' odd and m — m' even (NS 
representations) the different terms in K transform with the same phase up to ±1. Now 
let us compare the phase of Kx and K^. Since (^pp "^~^^\pp^^'^ '^'^^^ = —nm'p — for 
R representations (since mm' is even) K\ and transform with the same phase. For 
NS representations the phase of the terms in K\ can differ by a minus sign relative to 
the phase of the terms in Ky^. This is precisely the behavior expected from the characters 
because the Lq value of the different states in the representation differ by integer or half 
integer in the NS representations and they differ by an integer in the R representations. 
We conclude that under r ^ r + 1 the characters transform up to an overall phase as 
Xx^ ^ Xx^ cind Xx invariant. 

Now let us consider the behavior under r ^ — -. We use the Poisson resummation 

r 

formula 



a 

m 



Kx{t = ) = 2_^e ^pp'-' ^ \ ~ — 7 2^ ? 2pp' e 2pp' . (B.8) 



to write 



r ^ — ' V tpp' 

Consider for example the theory with {p = 2,p' = 8). It has two NS representations: 
the identity with A = 6, A = 10 and another representation with A = 2, A = 14. Using 



xn-r) = Xlln^oir) [Mr) - Kio(r)] = xfii,.=o(r) 
X^Hr) = xLlH=oir) [K2{r) - K^,{r)] = xfii..=o(r) 



q 32 _2^g 32 

n n 

q 32 -^q 32 



(B.9) 
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we easily find 



Xe I ~) — Xc=i,/i=ol ~) 

T T 



T T 



1 



j1 r 

I 32 



e 16 — e 16 



v^5/_i^ _ TVS 

X2 I -) — Xc=l,ft=0l ~^ 
T T 



lxfcW-;)Ef\/2^(9 



(16fc + l)"' (16fc + 7)'' 

32 — g 32 j 



2+ \/2(( 



i^2(--)-i^l4(--) 
T T 



(16fc + 3)'' 

r 32 



(16fc+5)-^ 
f 32 



g32 



g 16 — e 16 



2 + V2(g 



(16fc + l)^ 
32 



/ (16fc+3)^ 

+ V2- \/2(g^^ 



(16fc + 7)^ 
I 32 



(16fc+5)^ 
f 32 



(B.IO) 



Therefore, they can be expressed as hnear combinations of X&^{'^) ^^^d X2'^'^ (''")• 



Let us contrast this with a putative theory with [p = 2,p' = 6). It has two NS 
representations: the identity with A = 4, A = 8 and another representation with A = 0, 



A = 12. Now 



xr(r) = xfii,/.=o(r) [K,{t) - Ks{t)] = xfJi,H=oir) 
X^ir) = xf=l.=o(r) [Ko{r) - K^,{t)] = xfii,/.=o(r) 



(.en+iy 



E (6ri + 2)^ 
q . 



n 



(6ti + 3)^ 



n 



(B.ll) 
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and their modular transforms are 



-Xc=l,h=oi ~) 
T 



K,i--) - Ksi--) 
T r 



T 



^NS 



1, 



U=o(--)E^" 



j2. 

I 24 



27vin 27vin 

,6 — e 3 



(6fc + l)^ 

q 24 



(6fc+3)^ 
I 24 



-Xc=l,h=0\ ) 
T 



Ko{ — )-KM — ) 

T T 

n 



r 

T2i 



:b.i2) 



These are not hnear combinations of characters of the theory. Therefore, the theory with 
{p = 2, = 6) is not modular invariant. More generally, when p and p' are even we must 
assume that {p — p')/2 is odd jS^. 



Appendix C. Some properties of the ZS hierarchy 

In this Appendix we derive the form of the Zakharov-Shabat operators Rm, Hm which 
is relevant to solving string theory on a sphere. We also review how the mKdV hierarchy 
is related to the equation ( |6.5|) for the matrix model. 

We first recall that, in terms of Hm and Rm, the recursion relations are given by ( |2.13| ). 



It is interesting to solve the recursion relations (|2.13| ) on the sphere, i.e. by dropping the 
derivative terms and using the ansatz 



[m/2] 

(C.l) 

[(m-l)/2] ^ ^ 



m^2l+2^m-2l-l 



1=0 



When we insert this into the recursion relations ( p.l3| ) we can drop the term involving 
derivatives of H in the recursion relation for R^+i- Equating coefficients on both sides we 
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find the recursion relations for tlie coefficients: 

m+l 



bT^' =br - ' (C-2) 

(m - 2/)6J"+i = - (m - 2l)aY' + {m - 21 - 1)6^ . 

Note tliat tlie last two equations have to be compatible with each other. This implies the 
equation 

This relation can be stated as 

which is precisely the integrability condition for the second equation in (|2.13|) , when we 
drop derivative terms in Hn and Rm- Note also that once we set ag = 1 ( |C.2| ) implies that 
ttg^ = 1. By demanding that the right hand sides of the first two lines in (|C.2|) obey (|C.3| ) 
we find that 

By defining new variables p and (p through 

2 2,2 ^ . f /n c\ 

p = r + w , COS = — , sm (/? = -, (^-o ) 

p p 

we can see from ( |(J.5D that and i/^, can be written in terms of Legendre Polynomials 

• 2 

Hm=- P'"+MC0S<^P^(C0S(^) - P^+i(cOS<^)] = ^^P^ (cOS (/p) , 

Rm =p"'"^^ sin(^P^(cos(/?) , 

which is our main result. 

This relation can be derived more directly by looking at equation ( p.7| ). In the limit 
that r and uj are independent of x (i.e. commute with we have the equation 



0= I A ^= f- 'jt!li =i^J (C.8) 



2n ip + rJi + {uJ - QJs J 2n p"^ + \v\'^ / 4 



We see that 



Vl 1 ^ oo 

"T^ = 7T7 sincp 5 = — C~^~^ p'''^^ sin (zpP/ (cos 0?) 

N ICI ^(i-2cos(^f + f^)V2 |C|£^' ^ ^ u 

P (cos^J-^) ( °° 



^ (l-2cos.^;g)V. - ^ E C-V-[cos.PKcos.) - P.,.(cos.)] . 

(C.9) 
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For C < these equations imply (|C.7|) . Note that in this hmit the quantities O; in ( |2.10| ) 
are zero since they have no term without derivatives. 

Note that ( p.7| ) imphes that the resolvent, Tr[-p3^], is proportional to dx'Tr[JsO{x') 
with z = iC,. In the planar limit we can use (|C.9|) to compute it. We find that the ends of 
the cuts are aX z = iC, = ipe^'^'^ = ±r + iuj. 

Let us first consider even m and g = 0. The equations on the sphere are 

psinip[-amX + p'^Pmicostp)] = , 

(C.IO) 

p'^+^sin^ipP'^{cosip) =0 , 

with [—l)'^/'^am > 0. First note that sinip = is a trivial solution of the equations with 
vanishing r and vanishing free energy. Aside from this trivial solution, the second equation 
implies that P^{cosip) = 0. For a; > there is always a solution with cos 99 = 0: this is 
the symmetric solution with vanishing lo. The number of possible solutions actually grows 
with m. For m = 2 the solutions with cos f = 0,1 are the only ones. For m = 4 and 



X < we also find solutions with cos 99 = ± a/3/7: these are the broken symmetry solutions 
discussed in section 4 (see ( p.4|) ). For m = 6 and x < we again find two broken symmetry 



solutions with cos cp = ±Vl5 + 2v^/V33; for X > 0, in addition to the standard solution 
with cost/? = 0, there are also non-trivial solutions with coscp = 

general, the equation P^(cos(/?) = admits m — 2 non-trivial broken symmetry solutions; 
some of them are compatible with a: < and others with x > 0. For all the nontrivial 
solutions, p, r, and uj scale 

ii/m large Ixl. 

If m is odd then let us choose am > 0. For x > we look for solutions of the 
second equation in ( |C.10|) with Pm{cosLp) > (the solutions with x — > — x are obtained 



by cos If — cos f). The number of solutions for a given sign of x is (m — l)/2. After we 
include derivative terms in the string equations, the sphere solutions receive higher genus 
corrections. 

Consider, for instance, m = 3, which is the first non-topological "odd" critical point. 
Here the equations are Rs = 5xr/8; H3 = (see ( p.l4| )). We find the following solution as 
X — >• 00: 

r(x) = xi/3 - Ix-y^ - ^x-'> + 0(x-23/3) , 
9 729 

u{x) = -\x"' - - ^^-^ + 0{x-''/') , (C.ll) 

^ = ^ = -x2/3 - - ^x-i^/3 + 0(x-22/3) . 

(ix2 4 4 9 729 ^ ^ 
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Remarkably, for m large enough that there are multiple solutions, we find that the one loop 
partition function is independent of the choice of solution (except for the trivial solution 
r = where it vanishes). The results we find are consistent with the general formula 

Finally, we extend ( p.l4|) by presenting a few more terms in the Zakharov-Shabat 
hierarchy, generated via the recursion relations: 



R4 

R5 



3r" 



3r^ uj^ + ruj'' - 



51 ^ 
r r 



+ 12 ujr' ljj' + 3ruj'' 



5 r r 



2 ^if 



6uj'^ r" + 4ruj u" + 



—3 r uj 
2 

15 UJ 



+ 2r'^ uj^ — 2uJ r' + 2 r r' iv' + 4 r uj r" + r'^ uj 

2 



5r^uj^ + ruj^ 



25rujr'^ 25r'^r'uj' „„ 9 , , ^_ ,2 



25 r"^ ujr" 



3 II 



5 uj" 



2 2 

5/^(3) +5^r(^) +ru;(^) 



+ 30c<; r' uj' ^ UJ uj' — 



+ 10 r uj" + 10 r" uj" + 10 uj' + 



,2/2 



5r'^ Ibr'^uj'^ br'^uj^ 5r'^r''^ ,-2/2 / / br^uj'"^ 

--ns =-77^ ; 1 : buj r l\j r uj r uj -\- 



16 4 
5 r"^ r" 



2 4 

//2 

+ 10 r cc;^ r" + h 5 a; a;" - r' r^^^ + r r^^^ 

2 2 



(C.12) 



Appendix D. Double scaling limit of the complex matrix model 

In this appendix we derive ( |6.5| ) for the simplest case, /c = 1, where it reduces to (|6.7|). 



This is a review of the discussion in |3^], which is further expanded in ||32 



We start with the complex matrix model ( |6.1| ), which reduces to the integral (6.2) 



N 



'-V{yi) 



:d.i) 



^ ■ I'OO 

^ = n / dvivi^' 

We define orthogonal polynomials with respect to the measure 

//•oo ^ 
d^PnPm = J dyy'^e~~^ PnPm = Sn,mhn , Pn = 2/" + " 

D.l. q = 

Let us first set q = 0. We derive recursion relations in the usual way by writing 

yPn = Pn+l + SnPn + TnPn-l, = K/hn-1 , (D.3) 



(D.2) 
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and then writing 



N h 



n 5 



n 



where (m|n) = hnSmm and Xn = ^ and O are defined in terms of the values of the 
polynomials at zero 

They satisfy 

(D-6) 

r Q Q 1 = Q2 

'n'^'n'^'n — 1 '^'rt • 

The simplest potential, which arises for the k = 1 model, is V{y) = —y + and 
then we find that (|D.4|) are 

1 ~t~ Sri 7 • (■^•'^) 

We can eliminate Sn from this equation and substitute in ( p.6|) to end up with 



fn •^n 5 

+ 0^+1 = -(0^ + 1)0^ , (D.8) 
r Q Q 1 = Q2 

When the cut is far from y = 0, the polynomials are small at the origin; therefore, we find 
n = Vtn = and recover the standard hermitian matrix model equations. 

One can then combine these equations with ( p. 41 ) to obtain the string equations. In 
the simplest model this works as follows. Let us start by considering the planar limit. The 
equations ( p.6[ ) and ( p.7| ) (or equivalently ( p.8| )) have three solutions on the sphere: 

1. Q = il = 0, s = 1, r = X. 

2. O = i(-2 + VI + 12x), O = ^(1 - 12x + V1 + 12x), s = ^(1 + V1 + 12x) and 
r = ^ = ^{l + 6x + VI + 12 x). 

3. O = i(-2 - VI + 12x), O = - 12x - V1 + 12x), s = 1(1 - Vl + 12x) and 
r = ^ = ^(l + 6x — Vl + 12 x). At this point we need more physical input. It 
arises from ( p.5| ) and the interpretation of O in terms of Pn{0). First, it is clear that 
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n > 0. This rules out the third solution and allows the second solution only for x > |. 
Second, it is clear that we cannot have Pn(0) = for all n. Therefore we cannot take 
the first solution for all x. We conclude that for < x < ^ we should take the first 
solution and for | < x < 1 the second solution. In terms of the eigenvalues, for x < | 
they stay away from the origin so that = 0; but for | < x < 1 they reach the origin 
so that O 7^ 0. 

Now, consider the double scaling limit around the transition point x = 1/4, 

r = ^ - ew, X = ^ - e5, O = eO, O = -eTZ, N = 7e"^/^ (D.9) 

Equations (p.8| ) become 

TZ = u — z , 



2n - t^/^n' + -n" + 0{e^'^) = O + eO^ , 



(D.IO) 



where derivatives are evaluated with respect to z. Solving the second equation for O, and 
substituting into the last one, we find 

-nn'' -hn'f -An^{n + u) = Q . (d.ii) 

2 4 
Substituting 'u = 'u + 5/2we end up with 

^uR^ - ]-nn" + \{n'f = o , 

n = u-~z/2, (D-12) 
F" = 4u . 

D.2. q>0 

For g > 0, we start again with the left-hand sides in (p.4|) and integrate by parts. 
The boundary term vanishes, but there is an extra term where the derivative acts on the 
factor of A*^. We again call them O^^, 0„: 



1Q 
N ^ 



J y J y J y 

J y 

=/i;liPn(0)Pn-l(0) j d//^ + P^-l (0)P:(0)/lo//in-l 
=/i;llPn(0)Pn-l(0) j dtl^ + Pn{0)P;,_Mho/hn-l + 1 ■ 



(D.13) 
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Since the measure includes a factor of y'^, all these integrals converge. Note that in the 
formal limit — we recover the values for O and Q of ( p.5|) . Using (|D.13|) (and the two 
forms of O), we find the generalization of (p.6| ) 

The corrections due to q are of order therefore, they do not contribute in the 

planar limit. However, they have to be kept in the double scaling limit. Repeating the 
derivation of the differential equation ( p.ll| ) we find 



16:^^ + l67^2u + (7^')^ -2nn" = q^ 

TZ = u — z . 



(D.15) 



Substituting u = u + z/2 we end up with 



32n^u + {n'f -2nn" = q^ , 

n = u-z/2, (D.16) 
F" = 4u . 

Equations (|6.7| ), ( |6.4| ) may be brought to this form by defining 

u = 2^/^u , z = 2^/3^ , (D.17) 

from which it follows that TZ = 2^/^'R,. No rescaling of q is needed. This completes the 
derivation for k = 1. 

Note that (|D.13|) also implies that O > 0, and that it approaches zero as z —>■ +00. 
The equations ( p.6|) imply that in the scaling limit, where Sn is close to one, O is negative. 
Then ( p.9|) implies that TZ is positive and that 7^ ^ as z — > 00. It can also be seen from 
the second sphere solution above that 7t~ —5/2 for z ^ —00, which is indeed a property 
that we get from (|6.5[) once we impose that u —>■ as z ^ —00. 



D.3. The Miura transformation 

In this subsection we consider the ZS hierarchy for (7 = = w, so that f = r. In this 
case only the Fi in ( |2.7| ) with / even, / = 2k, are nonzero. These are the terms in the mKdV 
hierarchy. There is an interesting relation between these F2kif,g = 0) = R2k{f-,oJ = 0) 
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and the Gelfand-Dikii Polynomials Qk{u) [0. The functions can be defined in terms 
of matrix elements of the operator 



(rf, + /ji-CJ3) ' ^ ' - -'(^d^^fji-(:h){d:,-fji-(:h) 

1 (D.18) 

Hdx - fJi - CJs)— — 72 — — 

dl - f'Ji - V - CV4 

We can derive a similar relation where we write the numerator on the right. Taking the 
(a;| ■ ■ ■ matrix elements of these operators, extracting the piece proportional to 1, which 
should vanish, and the piece proportional to Ji, which should equal —J2^2kC~^'^~^j we 
find that 

F2k = 22^+i(c/, + f)Qkiu-) = -22^+i(d, - f)Qkiu+) (D.19) 

where 

These relations can be viewed as arising from supersymmetric quantum mechanics, where 
the operator O is the resolvent of the supercharge, and the operator d'^ + u is the Hamil- 
tonian. 

Let us introduce a derivative = d^ + f and assume that u in (|6.5|) can be written 
in terms of / through ( p.2Cl| ) as u = U-. Then it is possible to show that ( |6.5|) becomes 

q2 = Lnv+n - ^nvln + \{v+nf = -^nd^{v+n) + \{v+nf (D.21) 

Acting with on the second line in ( |6.5D we get 

V+n = ^(2/c + l)i2iF2i - 1 , i2k = (D.22) 

A;>0 

Note that in ( |6.5|) we are implicitly saying that to = — 4x. This translates into to = 
This X dependence of to is the source of the —1 in the right hand side of ( p.22| ). These 



relations were found in p9| , |33[] 



In conclusion suppose that we find a solution / that solves ( p.22|) with V^IZ = 2q 
with q >0. In other words, a solution of the equation 

2q + l = ^(2A; + l)t2fci^2fc (D.23) 

A;>0 
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Then we wiU also solve ( p.21| ). In ||2^ it was argued that smooth solutions exist for the 
mKdV equation with a zero left hand side. Since (|D.23|) follows from a lagrangian that 
is bounded belowS, if the highest nonzero t2k has the right sign, then it is easy to see 
that a smooth solution should exist that interpolates between / ~ x^/^'" for large x and 
/ ~ —{2q + l)/x for large negative x. Note that if 2^ + 1 is zero, then / —f is a 
symmetry of the problem. On the other hand if 2q + 1 is nonzero the minimum of the 
action at large x with positive / is the one with lowest energy. 

This solution is such that u obeys the physically relevant boundary conditions u 
^i/m ^ _^ g^j^^ w for X —> oo. We should also impose that 7^ > and that it 
goes to zero for large x, as we showed above. By integrating the equation Vj^R = 2q we 
find that 

- r f r r f 

n = 2qe ■'-0 / dye-^-o (D.24) 



For q > and with / obeying the above boundary conditions we see that TZ has the 
requisite properties. For q = it was proven in that a unique solution of exists. 

We can study the concrete example of /c = 1 in order to see how this works |]33,3l|. 
Defining u = U- in ( p.2Cl| ), we find 

2q + l = -\f" + -^f-xf (D.25) 

We see that this equation comes from an action bounded below and that the lowest mini- 
mum at large x has / ~ ^x^^'^. 

Note that if we had imposed instead Vj^IZ = —2\q\ we would also solve ( p.21|) , but 71 
would not obey the right conditions. 

This relation between ( [Ol ) and the equation ( p.23| ), can also be understood as follows. 
When we are dealing with a hermitian two-cut model with a symmetric potential, we can 
set g = and we can consider symmetric and anti-symmetric polynomials independently. 
Namely, the free energy can be expressed as a sum F = -f- F- where F± represents 
the contribution to the integral from even and odd polynomials. Then we have that 
Fj^ = 4" "F 4"' where / obeys the same equation for both cases. Of course, F" = /^/4 
is the equation for the free energy for the two cut model. In the hermitian model we can 



add a logarithmic potential of the form M log | A | [E9[ , which introduces a factor of | A | in 



We have checked up to A; = 3 that all the terms with derivatives in J i^2fe+i(/, </ = 0) are 
positive definite. We think it is true in general. 
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the measure that defines the orthogonal polynomials. The resulting equation for / is the 
same as in ( p.20| ) but with a constant proportional to M added to the right hand sideii 
(remember that we have set g = 0). Let us consider now the complex matrix model, and 
look at the first form of the integral in ( |6.2D . Then we can take the even polynomials 
P21 of the hermitian matrix model defined with the measure = A^+^^e"-"^" >. These 
polynomials are functions of and we can rewrite the Vandermonde determinant that 
appears in the first line of (|6.2| ) in terms of them. This will lead to an expression for 
the partition function of the complex matrix model in terms of the norms of the even 
polynomials of a two-cut model with a logarithmic potential with coefficient M = 1 + 2q. 
This last problem is precisely the one solved by the equation ( p.25| ) , with free energy given 
by F" = p/8 - r/A. 



Appendix E. The loop equation of the complex matrix model 

For M a complex rectangular x (A^ + q) matrix, with q positive, there are two 
closely related resolvents 

Let us derive their loop equations. We compute 

f dMdM^Tr ( -f-M-— 1 ^ ^-NT.viM^M) ^ q (^.2) 

J \dM M^M-zJ ^ ' 

which vanishes because it is a total derivative.^ This refiects the invariance of the integral 
under the change of variables 5M = Mj^^j^—^. It is straightforward to calculate ( [Li. 21 ) 

- ( Tr ^ ) + iN + q)( Tr — 7 ) - ( Tr — 7 Tr — 7 

\ M^M-z / ^ M^M-z/ \ M^M-z M^M 

^/^M^MV'(M^M)-zV'(z)\ 



M^M -z / ' ' \ M^M - 



+ q{Tr—-^ )-z(tt ^ 



M^M -z/ \ [M^M - zY 



{R{z)) 



Az 



(E.3) 



Normalizations in [|9| are different, fMinahan = fheref^- 
Compared to ( |6.1| ), we absorb I/7 into the definition of V{M^ M). 
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with f{z) = —jj^ (ti ^ MiA4^z ) ^ polynomial of the same degree as V'{z). 

In the large limit the term (^R[z)'^') factorizes. We do not neglect the term propor- 
tional to q = because we allow for the possibility that q is going to infinity with finite 
q. We derive the loop equation 

R{zf + {v\z) - ^) R{z) - ^ = (E.4) 

where we denoted the expectation value of the operator R by R. As a check, we can 
substitute i? = i? + | and derive the loop equation for R 

R(,f + {v'(.) + i\m.)-M^^^^^o (E.5) 

V z 4z 



i.e. exactly the same as ( [E.4|) but with R ^ R, q ^ —q, and the transformation of f{z) 
which follows from M^M ^ MM^. 
The solution of ( Fl ) is 



2R{z) = - (V'{z) + f ) ± y (V'iz) + f ) ' + ^ (E.6) 

We see that the parameter z in R{z) takes values in a two fold cover of the complex plane 
which is the Riemann surface 

y^=(^'W + f)' + ^, (E.7) 
or equivalent ly in terms oi y = yz 

f = {zV'{z) + qf + zf{z). (E.8) 
As z ^ oo on the upper sheet, we have R ^ This determines 



2R{P±{z)) = - (V'{z) + f ) ± y (v'iz) + f ) + ^ (E.9) 

where P±{z) denote the points on the upper and lower sheets which correspond to z. 
Finally, note that R has a pole on the lower sheet R ~ ~p"^7)- The other resolvent has a 
pole in the upper sheet R ~ r^ij ■ 
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Let us consider the simplest model with V{z) = + -2^/2). The polynomial 

y'^ in ( [li).(j| ) is of fourth order, describing a genus one surface. The polynomial f{z) is 
/ = —4z/'j + fo- If we are interested in a one-cut model we impose that the polynomial 
has a double zero. This determines the constant /q. The most general such polynomial is 

f = ^{z- a)2 [z'^ + 2z{a - 1) + (1 - 4a + 3a^ - 47 - 2-fq)] (E.IO) 

where a obeys the equation 

4a3 - Sa'* + -f^f + 0^(47 - 1 + 2-fq) = . (E.ll) 

We are interested in the double scaling limit, where one of the ends of the cut approaches 
2; ~ 0. Let us first set ^ = 0. From ([E.ll|) we see that a = is always a solution and leads 
to a cut extending from 1 — to 1 + This solution makes sense as long as 47 < 1, 
otherwise the cut would extend to negative values of z which are not allowed physically. 

Thus, for 47 > 1 the solution must have a 7^ 0. We write 47 = 1 — 5 and take the limit 
of small negative 6 and small a. From ( [E.llj ) with q = 0, we see that we get a solution 
with 4a = 5. For q nonzero we see that we need to scale q as Then we find that 

( PH.III) simplifies to 

4a^ + — - 5a2 = , (E.12) 
16 

and the second derivative of the free energy is proportional to the position z of the cut 
closest to the origin, which is given hj u = —4a + 5, in the scaling limit. After we scale 
-a ~ q'^^^v, and -6 ~ q^^H, then (|E.12| ) reduces to (|4.13|) . 



To summarize, for q = we have a phase transition between the regime with a = 
and a 7^ 0. In one phase y'^ has a double zero at z = 0, and in the other phase the double 
zero is at negative z while the cut reaches z = 0. For q ^ there is no such transition, 
and the cut does not reach z = 0. 

Let us now discuss the branes of OA theory. The FZZT branes with 77 = ±1 in the 
OA are the same as the FZZT branes with 77 = =i=l in the OB theory (at c = 0). We think 
that the resolvent ( [E.l|) corresponds to the FZZT brane with rj = — lEl and fx^ = —z. 
Indeed we see that the disk diagram has the expected one cut when expressed in terms 
of z ( |3.26| ). On the other hand, for negative /U, this brane will have an expectation value 
similar to the FZZT brane with 77 = 1 at positive \i ( |3.27| ), which in terms of z has a cut 
at the origin as expected. 
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Or the FZZT with = +1 in OB. 
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